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ABSTRACT 

We present rotation curves of 19 galaxies from THINGS, The H I Nearby Galaxy Survey. The high spatial and 
velocity resolution of THINGS make these the highest quality H I rotation curves available to date for a large 
sample of nearby galaxies, spanning a wide range of H I masses and luminosities. The high quality of the data 
allows us to derive the geometrical and dynamical parameters using H I data alone. We do not find any declining 
rotation curves unambiguously associated with a cut-off in the mass distribution out to the last measured point. 
The rotation curves are combined with 3.6 ijm data from SINGS {Spitzer Infrared Nearby Galaxies Survey) 
to construct mass models. Our best-fit, dynamical disk masses, derived from the rotation curves, are in good 
agreement with photometric disk masses derived from the 3.6 /im images in combination with stellar population 
synthesis arguments and two different assumptions for the stellar Initial Mass Function (IMF). We test the Cold 
Dark Matter-motivated cusp model, and the observationally motivated central density core model and find that 
(independent of IMF) for massive, disk-dominated galaxies, all halo models fit apparently equally well; for 
low-mass galaxies, however, a core-dominated halo is clearly preferred over a cuspy halo. The empirically 
derived densities of the dark matter halos of the late-type galaxies in our sample are half of what is predicted 
by CDM simulations, again independent of the assumed IMF. 

Subject headings: galaxies; fundamental parameters — galaxies: kinematics and dynamics — dark matter — 
galaxies: structure — galaxies: spiral — galaxies: dwarf 



1. INTRODUCTION 

Early observations of the rotation curves of spiral galax- 
ies clearly in dicated the pre s ence of "missing" matter in these 
objects (e.g.. lBosmalll978l: iRubin et al.lll978l: lBosmalll981l: 
see also lSofue & Rubinll2001h . This "missing" or "dark" mat- 
ter, as it is now more commonly called, has since become 
one of the pillars of modern cosmology. An elaborate frame- 
work, explaining and describing the properties of the Uni- 
verse has been built, involving a cosmological constant A, 
and a form of collisionless Cold Dark Matter (CDM). The 
ACDM paradigm provides a comprehensive description of the 
Universe at large (a s shown by, e.g., the WMAP results, see 
ISpergel et al.ll200'7h . but problems on smaller (galaxy) scales 
remain. 

A well-known problem is that galaxy-scale structures 
formed in cosmological ACDM simulations are too com- 
pact. One symptom of this is the well-known missing satel- 
lite problem (sir nulations produce an overabundance of dwarf 
galaxy halos, cf. lMoore et al.ll999al) . This seems to have been 
alleviated somewhat by the recent discovery of many faint 
dwarfs in our Local Group - some of which seem to be dom- 
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inated by dark matter - but the discrepancy is still significant 
(ISimon & Gehall2007h . 

The so-called cusp-core discrepancy (the predicted dark 
ma tter density profile s in galaxy halos are too steep; 
cf. Ide Blok et alJl200lT) is another manifestation of the (too) 
compact structures in CDM simulations. Galaxy rotation 
curves have in the last decade taken center stage in the debate 
about this problem. Going into some more detail, results from 
cosmological ACDM simulations suggest that the density pro- 
files of dark matter halos should be nearly universal, indepen- 
dent of the mass of the halo. This density profile has a charac- 
teristic steep mass-density slope in the inner parts which can 
be approximated with a power-law p ^ r" with a < -1; this 
is the so-called "cusp" (e.g., [Na varro, Frenk & WhitdlI996l 
I997|:lKlvpin et al.ll2001l;lHavashi et alj|2004l; iDiemand et alj 



20051) . 

Extensive observational determinations of this inner mass- 
density distribution seem, however, to indicate that the mass- 
density profiles of dark matter halos can be better described 
using an approximately constant-density inner "core" (p ~ r" 
with a ~ 0). T h is core has a typical s iz e of order a kpc (e.g.. 



Moord 
1997t 



U). 1 h is core has a typical s iz e or order a Kpc (e.g., 
T994I; Ide^k et all [19961 Ide Blok & McGaughl 



de Blok et al.1 120011: Ide Blok & Bosmal 120021; 



Marchesini et alj 120021; IWeldrake et al.l 120031; iGentile et alj 



2005t iKuzio de Narav et al.l 120061; iGentile et al.1 120071) . The 
discrepancy between these two outcomes has been much 
debated, as the presence of cusps is a fundamental feature of 
current cosmological simulations, confirmed many times by 
large increases in resolution and computing power that have 
been achieved over the years. 

A large part of this debate has focused on possible sys- 
tematic effects in the observations. Resolution of the data, 
accuracy of the positions of centers of galaxies, as well as 
presence of non-circular motions have featured prominently 
in these discussions (e.g., Swaters 1999; van den Bosch et alj 
2000; S waters, Madore .& Trewhella 2000; Simon et al. 200iF 
Rhee et aT] 120041) . Nevertheless, multiple, repeated. 
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independent long-slit measurements of rotation curves 
as well as analyses using two-dimensional optical ve- 
locity fields consistently yield results suggesting the 

presence of a dark m atter core in the inn er parts 

of d i sk galaxies (e.g., Marchesini et al. 2002; de Blo^ 
2004t Ide Blok & Bosmal i2002: Kuzio de Naray et al.. .200l 



Gentile et al.ll2007l: ISoano et al.ll2008h ^ 

H I synthesis observations of disk galaxies would, in princi- 
ple, be eminently suited to address this problem. The H I disks 
of galaxies have a higher filling factor than the often patchy 
Ha distributions used for optical rotation curves. Measuring 
accurate velocities and modeling them across the entire, two- 
dimensional, filled H I disks is therefore one of the strengths 
of radio synthesis observations. This means that possible ef- 
fects of non-circular motions, which may not always be rec- 
ognized in one-dimensional slit observations, can be analyzed 
within the context of the surrou nding rotating disk (lOh et alj 
I2OO8I; llYachternach et al.ll2008l) . Furthermore, a determina- 
tion of the position of the dynamical center can be made, en- 
abling a direct and independent comparison with the position 
of the photometric center (Trachternach et al. 2008). In short, 
the rotation curve, as well as the uncertainties that could affect 
it, can be constrained by the same observation. 

The one major disadvantage which hitherto hampered mak- 
ing the fullest possible use of H I data in these analyses, has 
been their modest angular resolution, especially when com- 
pared to optical observations: to directly address cosmologi- 
cally relevant questions regarding the distribution of dark mat- 
ter requires high spatial resolution. For example, unambigu- 
ously distinguishing between a core and a cusp, requires phys- 
ical resolutions better than ~ 1 kpc (de Blok 2004). Also, to 
determine the importance of non-circular motions within the 
disks of galaxies one needs to resolve the size scales asso- 
ciated with features that cause these motions, such as bars, 
spiral arms and oval distortions. 

The majority of published H I synthesis observations have 
at best beam sizes > 15", meaning they can reach the desired 
sub-kpc resolution only out to distances of ~ 4 Mpc. This, 
of course, severely limits the number of galaxies available 
for study. THINGS (The H I Nearby Galaxy Survey) has ad- 
dressed this by observing in H I a sample of 34 nearby galax- 
ies at the desired sub-kpc resolution. To achieve this, it used 
the NRAO*^ Very Large Array (VLA) in its B, C and D con- 
figurations. The total observing time (including data retrieved 
from the archive) was ~ 500 hours. The maximum angular 
resolution that can be achieved with these data is ~ 6" for ro- 
bust weighting and ~ 12" for natural weighting. The velocity 
resolution is 5.2 km s~' or better With the sample galaxies 
at distances between 2 and 15 Mpc, the typical linear resolu- 
tion varies between ~ 100 and ~ 500 pc. For further discus- 
sion on the sample, aims and scientific goals of THINGS, as 
well as the technic al details of the observations, we refer to 
IWalteret al.l(l2008h . 

The current paper focuses on the rotation curves of 
THINGS galaxies. We use a sub-sample of the THINGS 
galaxies (discussed in Sec. 2) and present their velocity fields, 
tilted-ring models and the rotation curves. Additionally, 

' Note that the paper by Swaters et al. (2003), which is often quoted as an 
observational counter-example, states that their results "cannot rule out steep 
slopes", but also that "halos with constant density cores provide somewhat 
better tits". 

* The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 



we combine our Hi data with 3.6 fim dat a from SINGS, 
the Sp itzer Infrared Nearby Galaxies Survey (IKennicutt et alJ 
l2003h and derive well-constrained stellar mass-to-light ratios 
for the stellar disks. We use these to investigate the mass 
distributions in the sample galaxies using core- and cusp- 
dominated dark matter halo models. A number of galaxies 
presented here have been observed in the past, in some cases 
as part of the pioneering work on H I rotation curves in the 
1970s and 1980s'^, but this is the first time that such a large 
sample of galaxies is analyzed and modeled in a uniform way 
at sub-kpc resolution. 

In Sect. 2 we summarize the selection of our sample. Sec- 
tion 3 describes the creation of the velocity fields as well as 
the derivation of the rotation curves. Section 4 presents the 
rotation curves and velocity fields of the individual galaxies 
in more detail. In Sect. 5 the mass models are presented, and 
in Sect. 6 these models are discussed in more detail for in- 
dividual galaxies. The results are discussed in Sect. 7 and 
summarized in Sect. 8. Technical aspects of the analysis are 
presented in Appendix A. 

2. THE SAMPLE 

Our goal in this paper is to derive rotation curves of 
THINGS galaxies which depend on a minimum amount of 
external information. That is, we aim to derive the geomet- 
rical and kinematical parameters of the galaxies (inclination, 
position angle, rotation velocity, etc.) as much as possible 
from the THINGS Hi data alone. This limits our analy- 
sis to rotationally dominated galaxies with favorable orienta- 
tions. Although using additional information (such as derived 
from optical data) would, in principle, increase the number 
of THINGS galaxies for which a rotation curve can be deter- 
mined, the unquantifiable systematic effects associated with 
this introduce their own difficulties, as discussed below. 

There are practical upper and lower limits to the incli- 
nations of galaxies for which one can reliably derive rota- 
tion curves using standard methods. At very high inclina- 
tions (/ > 80°) the line of sight crosses a large range of 
projected velocities and it becomes difficult to quantify the 
typical rotation velocity (alth ough not impossible, see e.g., 
iKregel & van der Kruitl200l . Such high-inclination galaxies 
are, however, not present in the THINGS sample. 

At low inclinations , the practical limit lies at ~ 40° 
(lBegemanlll987l [19891) . At lower inclinations the velocity 
dispersion of the gas starts to become an important contami- 
nant. As rotation takes place in a two-dimensional disk, the 
line-of-sight (projected) component of this rotation velocity 
decreases in amplitude with decreasing inclination, whereas 
the contribution due to random motions, which take place in 
three dimensions, remains the same. This leads to unaccept- 
ably large uncertainties in the deprojected rotational velocity. 

Independent of inclination, we also require our sample 
galaxies to be dominated by rotation. This is not the case 
for a small number o f THINGS galaxies. These are M81dwA 
dSargent et al.ll983h and M81dwB as well as NGC 15 69 (dy- 
namics disturbed by a starburst; IStil & Israel l2002h. N GC 
3077 (heavily interacting wit h M81: IWalter et al.ll2002h and 
NGC 4449 (tidally disturbed; [Hunter et alj|19'99r We do not 
consider these galaxies further. 

To quantify our sample selection further, we consider in- 
dicative H I inclinations derived from ellipse fits to the outer 



' See, e.g.. Fig. 1 in lBosmj 098 Ih for a proto-THINGS depiction of the 
velocity fields of 22 nearby spiral galaxies. Half of these are also part of 
THINGS. 
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Hi disks, in addit ion to optical inclinations from LEDA 
(iPaturel etalj|2003h . These values are listed in Tabled In 
determining the H i inclinations, using the maps presented in 
IWalteretal.l(l2008h . care was taken that these fits were not af- 
fected by the tails, clouds, and other features often seen in the 
outer parts of the H I disks. The H I inclinations listed here are 
deemed to be as reliable as their optical counterparts, and in 
some cases even prefer able (cf. the values for DDO 154 and 
the associated images in lWalter et al.l200"8l) . We stress that de- 
termining precise values for these inclinations is not the issue 
here, as we merely use them to clarify our sample selection. 
They are not used for any subsequent dynamical analyses. 

A first analysis of the velocity fields presented in 
IWalter et al.1 (|2008|) showed that THINGS galaxies with in- 
dicative H I inclinations higher than 40° are all suitable for 
rotation curve analysis. These galaxies are all discussed in 
the current paper Galaxies with H I inclinations less than 30° 
are unsuitable and will not be discussed any further here. Of 
the galaxies with Hi inclinations between 30° and 40°, only 
NGC 6946 proved suitable. The reason why for this galaxy a 
rotation curve can be derived, whereas this was not possible 
for, e.g., NGC 4214 (which has a comparable indicative Hi 
inclination), is the large ratio of galaxy size to beam size. The 
resulting larger number of resolution elements in the NGC 
6946 velocity field help constrain the kinematics better than 
for the other galaxies in that (H l) inclination range. 

A look at Table [T] shows that there are a few cases where 
the H I and optical inclinations differ substantially. These can, 
however, be readily explained. DDO 154 (/hi > 40°, /leda < 
40°) has a well-defined and regular H I disk, but a low surface 
brightness ill-defined optical component. The H I inclination 
is thus far more reliable. M81 dwB (/hi < 40°, /leda > 40°) 
has a high-surface brightness central bar-like structure which 
affects the optically derived inclination. NGC 3351 (/hi < 
40°, /leda > 40°) has a fairly face-on and regular H I disk, 
but optically is dominated by a very strong central bar. NGC 
4736 (/hi > 40°, /leda < 40°) is the only galaxy where the 
outer H I disk seems to have a genuinely higher inclination 
than the central optical component'". 

One might argue that our choice to only consider galax- 
ies for which the geometrical parameters can be kinematically 
determined is needlessly strict, and that rotation curves could 
in principle be derived for the lower inclination galaxies by 
using additional information, such as the axis ratios of the op- 
tical or H I disks. However, in these cases, the best one can do 
is assume a constant inclination. This will then ignore intrin- 
sic inclination changes (due to, e.g., the presence of warps), 
which can strongly affect the derived rotation velocities: an 
estimate which is only a few degrees off at these low inclina- 
tions readily leads to spurious rotation velocities. To illustrate 
this, at an inclination of ^ 20°, a 5° uncertainty, equivalent to 
a 3 percent uncertainty in the axial ratio, leads to a 25 percent 
uncertainty in the deprojected rotation velocity. Furthermore, 
as discussed earlier, at low inclinations random or streaming 
motions can be as important as the projected rotation velocity 
and these will therefore severely affect the resulting rotation 
curve. 

In the following, we will limit ourselves to the subsam- 
ple of galaxies indicated in Table [T] i.e., all high-inclination 
THINGS galaxies (/ > 40°) that are dominated by rotation. 

This is consistent with the results we will derive later in this paper NGC 
4736 shows a well defined inclination trend, varying from / ~ 30° in the inner 
parts to I ~ 50° in the outer parts. 



and, in addition, NGC 6946, leading to a total sample size of 
19 galaxies. In our analysis we will use the natural-weighted 
THINGS data sets. The robust weighting does give a higher 
angular resolution, but at a decreased column density sensitiv- 
ity. We found that the small loss in resolution due to the use 
of the natural-weighted data is more than compensated for by 
the increased column density sensitivity and, hence, increased 
usable area of the velocity fields. 

We did apply neither residual-scaling corrections, nor pri - 
mary beam corrections (as described in IWalter et alj|2008l) . 
The reason for not using these corrections is that for the con- 
struction of our velocity fields (described in the next sec- 
tion), we need a constant noise level throughout each galaxy's 
data cube. Primary beam corrections introduce a noise and 
flux scaling dependent on distance to the pointing center 
Residual-scaling, in essence, artificially decreases the noise 
in the data cubes in order to produce correct flux densities. 
All velocities mentioned in this paper are heliocentric using 
the optical velocity definition. 

3. VELOCITY FIELDS AND ROTATION CURVES 
3.1. Velocity field types 

Hi rotation curves are most often derived from velocity 
fields. A velocity field aims to give a compact and accurate 
"short-hand" description of the dynamics of a galaxy by as- 
signing a "typical" velocity to every spatial position. That 
is, for every position one uses the velocity that most accu- 
rately represents the circular motion of the bulk of the quies- 
cent component of the gas as it moves around the center of 
the galaxy. In an ideal case (quiescent, purely rotating disk, 
infinite resolution, intermediate inclination) this velocity cor- 
responds to the velocity associated with the peak emission in 
the velocity profile. 

In practice, however, the observed velocity profiles are af- 
fected by several systematic effects. These can be divided into 
effects inherent to the observations themselves, and those due 
to physical processes within the galaxy. 

The first category, instrumental and observational effects, 
is especially important in lower resolution observations. Here 
the effect of the relatively large beam (compared to the size 
of the galaxy) causes asymmetric velocity profiles, with the 
long tail always pointing towards the systemic velocity. This 
is known as "beam smearing". Generally, one attempts to cor- 
rect for this by choosing a typical velocity close to the max- 
imum velocity (i.e. furthest away from the systemic veloc- 
ity) found in the velocity profile (e.g. Sancisi & Allen 1979). 
When left uncorrected, these resolution-related effects (in- 
cluding finite velocity resolution) always lower the apparent 
velocity. 

The second category is due to physical processes within the 
galaxy which cause non-circular motions. These will show up 
as deviations both towards and away from the systemic veloc- 
ity. For our galaxie s, this is discuss ed ext ensively in accom- 
panyin g papers by lOh et alj (l2008h and iTrachternach et all 

mm . 

To minimize the impact these effects can have, different 
methods to construct velocity fields have been used in the 
past (discussed below). The THINGS observations have high 
enough spatial and velocity resolutions t hat t he resolution ef- 
fects are no longer a problem (cf. Sect. 13.41 ): they also give 
us a detailed enough view of the small-scale movements of 
the H I within the galaxies to enable us to quantify the non- 
circular motions. Nevertheless, each method determines the 
"typical" velocity in a sUghtly different way, and it is worth 
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TABLE 1 

Indicative global inclinations of the THINGS galaxies 



Name /hi (°) 'LEDA (°) Name /hi (°) /leda (°) 



NCjL oio 


15 


21 




5U 


58 


NCjC 1569 




55 




d5 


71 


NGC 2403 


55 


61 


NGC 2841 


69 


61 


NGC 2903 


66 


64 


NGC 2976 


54 


58 


NGC 3031 


59 


58 


NGC 3077 




35 


NGC 3184 


29 


12 


NGC 3198 


72 


74 


NGC 3351 


39 


53 


NGC 3521 


69 


58 


NGC 3621 


62 


66 


NGC 3627 


61 


63 


NGC 4214 


38 


39 



NCjC 4449 




54 


JMjC 4 /jo 


44 


29 


NGC 4826 


64 


58 




5 1 


53 


NGC 5194 


30 


32 


NGC 5236 


31 


21 


NGC 5457 


30 


12 


NGC 6946 


35 


17 


NGC 7331 


77 


67 


NGC 7793 


43 


56 


IC 2574 


51 


67 


DDO 53 


33 


46 


DDO 154 


70 


29 


Hoi 


27 


17 


Ho II 


31 


45 


M81 dwA 


27 




M81 dwB 


28 


63 



Note. — Galaxies that have rotation curves presented in this paper 
are printed in bold. For some galaxies a representative inclination could 
not be determined. All inclinations listed here are indicative and are not 
used in the dynamical analysis presented in this paper 



exploring these to quantify their impact on the final results. 

In principle, our data are of high enough quality that (bar- 
ring large non-circular motions) the velocity of the peak emis- 
sion in a profile is representative of the total rotational motion. 
The most direct way to create a velocity field is therefore 
to simply determi ne the velocity at which this peak occurs. 
lUsero et akl ( l2008h derive these velocity fields in their study 
of the distribution of brightness temperatures within galaxy 
disks. For our purposes, these velocity fields are less suited, 
due to the decreased signal-to-noise in the outer parts of the 
galaxies and their smaller usable areas. As mentioned, in the 
following we briefly discuss alternative methods most com- 
monly used to construct velocity fields, list their pros and 
cons, and present our favored method. 

A practical comparison of the outcomes of the various 
methods is shown in Figs. [T]|4] and discussed in more detail 
below as well as in Sect. 3.2 and 3.3. 

First-moment or intensity-weighted mean velocity: One 
of the most well-known and widely-used methods involves 
taking the first moment or intensity-weighted mean along the 
velocity axis in the data cube. The procedure is straightfor- 
ward, but has the disadvantage that, if the velocity distribu- 
tions are not symmetric, the "typical" velocities are biased 
towards the longest tail of the velocity profile. Also, to avoid 
noise at outlying velocity values affecting the mean, one needs 
to carefully identify the regions with significant emission in 
the data cube and use only those areas when constructing the 
velocity field. This method is the one most commonly imple- 
mented in astronomical software. 

Peak velocity fields: Using the velocity where the H i pro- 
file flux peaks has the great advantage that it is totally inde- 
pendent from assumptions on profile shapes. A disadvantage 
is that as the profile peak value decreases the velocity values 
become more and more affected by noise. The peak velocity 
field is therefore best used in regions of high signal-to-noise. 

Gaussian profiles: An alternative is to fit single Gaussian 
functions to the velocity profiles. The velocity of the peak of 
the Gaussian profile is then taken to be the velocity field value. 
This method is less sensitive to noise and (moderate) asymme- 
tries in the profiles than the first-moment method, but can still 



suffer from biased velocities if the asymmetries are strong or 
if multiple components are present. The single Gaussian fit 
method is best used in cases where the instrumental resolu- 
tion is comparable to the typical width of the profiles. Many 
of the "classical" Hi rotation curve studies from the 1970s 
and 1980s fall in this category, and used this method to derive 
rotation curves. 

Multiple Gaussian profiles: A logical extension is to fit 
multiple Gaussian functions to the velocity profiles. In many 
cases, however, the narrow width of the H I profiles com- 
bined with limited velocity resolution yields only a relatively 
small number of independent profile data points with suffi- 
ciently high signal-to-noise. Simultaneous fits with two or 
more Gaussian functions can in these case become sensitive 
to the choice of initial values, or have more free parameters 
than data points. A careful fine-tuning is in any case needed 
to provide optimum fits. An additional uncertainty is how 
to decide which of the components represents the true rota- 
tion velocity, and which one characterizes non-circular, ran- 
dom motions (induced by, e.g., star formation or spiral arms). 
As stated, these uncertainties apply to simultaneous fitting of 
multiple Gaussian functions. An alternative method is to op- 
timally fit a single Gaussian function to the profile, subtract it, 
and subsequently make a second Gaussian fit to the residual 
profile. A robu st application of this method is presented in 
lOh et al.l ( l2008h . This method is remarkably effective in sepa- 
rating random motions from the average rotation velocity (or 
'bulk' velocity, to use the terminology from[Oh et al. 2008). 

Hermite I13 polynomials: A robust way to minimize the 
number of components is to include in the fitting function a 
prescription of the typical asymmetry encountered. A con- 
venient function to use for this purpose is a Gauss-Hermite 
polynomial that includes an /13 (skewness) term (see e.g., 
Ivan der Marel & Franxlll993b . The Hermite method has al- 
ready been applied s uccessfully to derive rotat ion curves of 
early-type galaxies in lNoordermeer et"an (l2007h . and we will 
use it as our primary method to derive velocity fields. We do 
not include an additional term in the fit, which would pro- 
vide a measure of the kurtosis of the profile (i.e., whether the 
profile is "fatter" or more "pinched" than a Gaussian func- 
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tion). Attempts to include an /14 term siiow that, for the pur- 
poses of determining a velocity field, this leads to diminish- 
ing returns; the fits become less stable and need increasingly 
fine-tuned initial estimates. For high signal-to-noise profiles, 
the velocities derived using the hermit profile will be close 
to those as derived from the peak-velocity fields. One could 
therefore argue that the latter are therefore preferable, as they 
are model-independent. However, in low signal-to-noise re- 
gions the Hermite velocity field results are more stable as 
more data points (from neighboring channels) are used to 
minimize the effects of the noise. We therefore prefer to use 
the Hermite velocity fields. 

We illustrate the various methods for two representative 
galaxies. The first is NGC 2403 (Sect. 3.2, Fig. [B, which 
is dominated by rotation, and where most of the asymmetric 
profiles are associated with streaming motions along the spi- 
ral arms. The second example is the dwarf galaxy IC 2574 
(Sect. 3.3, Fig.|2]i where a large fraction of the profiles is af- 
fected by asymmetries and random motions. 

3.2. NGC 2403: a galaxy dominated by regular rotation 

Figure [U shows the velocity fields for NGC 2403, derived 
using the procedures described above, together with differ- 
ence velocity fields with respect to the Hermite velocity field. 
The absolute differences between the Hermite velocity field 
and the Gaussian and first-moment velocity fields are larger 
than 10 km s~' for only a small fraction of the area. For the 
peak velocity field the difference with the Hermite field is con- 
sistent with pure noise (mostly caused by the lower signal-to- 
noise of the peak velocity field). 

Figure [3] shows an example of an asymmetric profile in 
NGC 2403, together with the values for the typical velocity, 
again derived using the methods described above. Shown in 
Fig. Hare the rotation curves derived from the respective ve- 
locity fields, using otherwise identical assumptions and meth- 
ods (cf. Sect. 3.5). Any difference is thus entirely due to the 
manner in which the velocity field was derived. The curve 
based on the Hermite velocity field systematically shows the 
highest rotation velocities and agrees best with the peak ve- 
locity field curve. The curves based on the Gaussian and 
first-moment fields are progressively more affected: the lat- 
ter curve systematically underestimates the rotation velocity 
by ~ 4 km s"' over the inner half of the disk. The bulk veloc- 
ity field curve is also consistent with the Hermite field rotation 
curve, but with a larger scatter. 

3.3. IC 2574: a galaxy with substantial random motions 

A second example is IC 2574, where the velocity field 
is clearly affected b y non-circular, random motions (cf. 
I Walter & Brinksll999h . These are clearly visible as "kinks" in 
the iso-velocity contours in all of the Hermite, Gaussian, first- 
moment and peak velocity fields shown in Fig. |2l It is clear 
that before a rotation curve can be derived some correction for 
the non-circular motions must fir st be m ade. This process is 
described in detail for IC 2574 in lOh et a l. (2008). The effect 
of this correction is illustrated in Fig. [3] The resulting "bulk" 
velocity field, showing the rotational part of the kinematics, 
is also shown in Fig. ID Though noisier than the other veloc- 
ity fields, the kinks have largely disappeared, and the velocity 
field is consistent with simple solid-body rotation. The result- 
ing bulk rotation curve also differs significantly from those 
derived from the uncorrected velocity field as shown in the 
bottom panel of Fig. |4] 



Further analysis shows that in our sample only NGC 2366 
and IC 2574 are significantly affected by these random mo- 
tions and they are therefore the only two galaxies where the 
bulk velocit y field differs s ignificantly from the Hermite ve- 
locity field (lOh et al.ll2008h . For these two galaxies we will 
in the rest of this paper use the latter. For the other galax- 
ies the differences are negligible and we will use the Hermite 
velocity fields. 

3.4. Resolution Effects 

The THINGS data are extremely well resolved. For exam- 
ple, the ratio of galaxy size to beam size for NGC 2403 is 
^ 350:1. Nevertheless, some galaxies, including NGC 2403, 
show steep increases in velocity in their inner parts, and de- 
spite the high resolution, these small radii could potentially 
still be affected by beam smearing. We quantify this effect 
by constructing two model data c ubes using the Groningen 
Image Processing SYstem (GIPSY: Ivan der Hulst et alJI 19921) 
task GALMOD. This task takes an input rotation curve, ra- 
dial H I distribution, inclination and position angle, as well as 
H I scale height and velocity dispersion and creates a model 
data cube at arbitrarily high resolution by distributing a large 
number of H I "clouds" through this data cube using the input 
parameters as probability distributions. This cube can then be 
"observed", e.g., by smoothing it to a lower resolution, and be 
used to construct velocity fields, derive rotation curves, etc. 

To make the model galaxies as realistic as possible, we 
adopt the H I distribution and average inclination and posi- 
tion angle of NGC 2403. We assume a velocity dispersion 
of 8 kms"', and a vertical (Gaussian) scale height of 100 pc. 
None of these assumptions critically affect the results. For 
the input rotation curves we adopt two extreme versions of 
the observed NGC 2403 rotation curve. For the first model 
we assume a steep linear rise to 130 kms"' (the maximum 
rotation velocity of NGC 2403) within 1 kpc (65"), and a flat 
130 kms"' rotation curve outside 1 kpc. The second model 
is identical, except that it rises to its maximum rotation ve- 
locity within 0.5 kpc (33") (compare these with the observed 
rotation curve shown in Fig.|2l and indicated in Fig.|5]). 

The resulting model cubes were smoothed to a resolution 
of 8" (the natu ral- weighted beam size for the NGC 2403 ob- 
servations; cf. I Walter et al.ll2b08h and a channel spacing of 
5.2 kms"' was adopted. The resulting "observed" cubes were 
then used to construct Hermite velocity fields and rotation 
curves were derived (cf. Sect. l3.5l and l3.6l ). 

We compare the derived curves with the respective input 
models in Fig. |5] Though we did calculate the cube for the 
entire radial extent of our pseudo-NGC 2403, we only show 
the inner portions of the rotation curves in Fig.|5]to focus on 
the rising part of these curves. It is clear that beam smear- 
ing is not a serious problem, with most of the curves differing 
by less than 1 kms"' from the input model. The only signif- 
icant difference occurs around R = 65" for the 1 kpc model 
and around R = 33" for the 0.5 kpc model, but this is entirely 
due to the unrealistically sharp break in our input curves. The 
innermost point of both curves is also somewhat affected, but 
only at the level of ^ 3 kms"' for the 1 kpc model and at 
^ 6 kms"' for the 0.5 kpc model. In summary, this demon- 
strates that resolution effects such as beam smearing have no 
significant effect on the THINGS data. 

3.5. Deriving the velocity fields 

As mentioned before, we use data sets without primary- 
beam and residual-scaling corrections as we need the original 
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Fig. 1. — Different types of velocity fields derived f(f SPjP^4()3. In the left column are shown, ^i^fflPft^ to bottom, the Hermite velocity [Her], the single 
Gaussian [Gau], the intensity-weighted first-moment (Moml), the bulk [Bulk] and the peak amphtude [T(peak)]) velocity fields. The velocity contours run from 
10 kms"' (top-right; light grayscales) to 250 kms"' (bottom-left; dark grayscales) in steps of 30 kms"' . The minor axis contour is at 130 kms"' . The right 
column shows the difference velocity fields with respect to the Hermite velocity field, as indicated in the figure. Grayscales run from —15 (black) to -1-15 (white) 
kms"', in steps of 5 kms"'. The white contour shows the —10 kms"' level, the black contour the -I-IO kms"' level. 



THINGS rotation curves 



7 



noise properties of the data to perform the profile fitting. For 
each galaxy, Hermite polynomials were fitted to all veloc- 
ity profiles with the GIPSY task XGAUFIT. 

To ensure high-quality velocity fields, we only retain pro- 
files with fitted intensity maxima higher than 3(7chan (where 
Cchan is the average noise in the line-free velocity channels 
in the relevant cube) and where the equivalent Gaussian dis- 
persion of the fit is larger than the channel separation (see 
IWalter etalJl2008l for the noise values and channel separa- 
tions). To objectively remove a small number of noise pix- 
els admitted by these filter criteria, we use the integrated H I 
column density map as an additional mask and only retain the 
fits at positions where the total flux value in the integrated H I 
map is higher than 3(Tn, where is the noise in the inte- 
grated H I map. For independent channels (as is the case with 
the THINGS cubes) aff is defined as = VNachnn, where 

is the number of channels with signal contributing to each 
pixel". These various steps result in objectively defined Her- 
mite velocity fields that form the basis of our analysis. 

3.6. Deriving the rotation cun'es 

The high spatial and velocity resolutions of the THINGS 
observations, as well as the favorable inclinations of the sam- 
ple galaxies, make these objects ideal targets for a tilted-ring 
rotation curve analysis. In a standard tilted-ring analysis a 
galaxy is described using a set of concentric rings, each with 
their own inclination /, position angle PA, systemic velocity 
Vsys, center position (xcjo) and rotation velocity Vc- Assum- 
ing that the gas moves in purely circular orbits, one can then 
describe the line-of-sight velocity for any position {x,y) on a 
ring with radius R as 



Vix,y) = V,y, + Vc(R)sm(i)cos{0). 



(1) 



In this equation, 9 is the position angle with respect to the 
receding major axis measured in the plane of the galaxy. This 
quantity is related to the position angle PA in the plane of the 
sky by 



cos(6') = 



sin(6') : 



-(x - xo) sin(PA) + (y - yo) cos(PA) 
R 

-(x - xq) cos(PA) - (j - yo) sin(PA) 
Rcos{i) 



(2a) 
(2b) 



PA is defined as the angle measured counter-clockwise be- 
tween the north direction on the sky and the major axis of 
the receding half of the galaxy. For each ring the parameters 
are varied using a least-squares algorithm until an optimum fit 
with the velocity field is achieved. 

We use the GIPSY task ROTCUR to make the tilted-ring fits 
and derive the rotation curves. As positions closer to the major 
axis carry more rotational information than positions near the 
minor axis, we use a | cos(0)| weighting for all fits. In defining 
the radii and widths of the rings we sample the rotation curve 
at a rate of two points per synthesized beam width. 

In this section we give a general description of the pro- 
cedure we used. Detailed, and more technical, descriptions 
of the derivation of individual rotation curves are given in 
the Appendix. For each galaxy we firstly fix the position 
of the dynamical center to that determined in the analysis of 

' ' Note that in creating an integrated H I map one adds emission across a 
number of cliannels, rattier tlian calculating tlie average. Provided the chan- 
nels are independent, i.e., the noise is not correlated between channels, the 
noise thus increases as the square root of the number of channels added. 



iTrachternach et all ( l2008h who used the same data. For com- 
pleteness, the resulting coordinates of the dynamical centers 
are given in Table |2l After fixing the center position, an ad- 
ditional ROTCUR run is used to determine Vsys (with PA and 
/ still left as free parameters). In general the systemic veloc- 
ity has very little uncertainty associated with it. The adopted 
values are given in Table |2l After fixing the central position 
and systemic velocity, ROTCUR is run once more with PA and 
/ as free parameters. The trends of / and (especially) PA with 
radius are generally well-defined, varying smoothly with ra- 
dius and exhibiting little scatter. Of the two fitting param- 
eters, the PA is very stable, and any small perturbations in 
its value result in only second-order changes in the rotation 
curve. In some cases, small-scale fluctuations in the incli- 
nation are seen. These are mostly caused by effects such as 
streaming motions along spiral arms which the fitting program 
tries to "compensate" for by changing /. 

To prevent these small spurious changes and to retrieve the 
underlying "bulk" rotation, we describe the PA and ; distri- 
butions by slightly smoothed radial distributions that ignore 
these small-scale "wiggles". A simple box-car smoothing 
with a kernel width of 3 or 5 data points was used in most 
cases. This only affects the point-to-point scatter of the / and 
PA values for the few galaxies where these wiggles are rele- 
vant, and does not affect the resolution of their radial distri- 
butions. With PA and / fixed using these radial distributions, 
we derive the final rotation curve. We also construct model 
velocity fields based on the tilted-ring models, as well as cor- 
responding residual velocity fields containing the non-circular 
component. 

The inclination wiggles described above, as caused by 
streaming motions, affect only a small part of the area of the 
velocity field for a small number of galaxies and do not af- 
fect the shape of the rotation curves in any critical way. Our 
treatment of the ; and PA distributions assigns the dominant 
non-circular motions to the residual velocity field. An explicit 
comparison between the ROTCUR residual velocity fields and 
the results from a full rigorous treatm e nt of the non-circular 
motions is made in ITrachternach et"al1 (l2008h and shows that 
the above procedure is justified. Care was also taken to 
double-check that the solutions produced using our PA and 
; distributions were self-consistent (i.e., we checked whether 
running ROTCUR again with only PA and / free, and Vq fixed, 
yielded PA and ; distributions consistent with our input val- 
ues). 

The method just described assumes galaxies are az- 
imuthally symmetric. However, asymmetries in the disk, or 
lopsidedness, can introduce intrinsic differences. To quantify 
the uncertainties introduced by this, we derive separate fits to 
the rotation curves of the approaching and receding sides of 
the galaxies. 

A final issue is the definition of the final uncertainties in 
the rotation curves. Unfortunately, there is no consensus in 
the literature on how to quantify these. Sometimes the for- 
mal Icr x^-fit uncertainty in Vc is used. However, by their 
nature, these uncertainties are much smaller than the disper- 
sion of individual velocity values found along a tilted ring. 
This dispersion is, therefore, also commonly used as an alter- 
native definition, as it is more representative of the physical 
uncertainties than the formal fit error For the NGC 3198 ro- 
tation curve (discussed in Sect. 4.8), for example, the average 
1(7 fit error is 17 times smaller than the average disper- 
sion along the rings. Using only formal x^-based error bars 
thus severely underestimates the true "physical" uncertainty. 
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Fig. 2. — Different types of velocity fields derived ftffl5?T4. In the left column are shown, ffonPE^!^o bottom, the Hermite velocity (Her), the single 
Gaussian (Gau), the intensity-weighted first-moment (Moml), the bulk (Bulk) and the peak amplitude [T(peak)]) velocity fields. The velocity contours run from 
—5 kms"' (bottom-right; light grayscales) to +100 kms"' (top-left; dark grayscales) in steps of 15 kms"'. The right column shows the difference velocity fields 
with respect to the Hermite velocity field, as indicated in the figure. Gray-scales run from -15 (black) to +15 (white) kms~\ in steps of 5 kms~'. The white 
contour shows the —10 kms"' level, the black contour the +10 kms"' level. 
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Fig. 3.— Left: Observed profile of the H I emission in NGC 2403 at position (a, 5)(2000.0) = 07''37"'36.7',65°36'46.3" (thick line). We adopt the notation 
given in Fig.[T] The five vertical lines indicate the derived typical velocities using, from top to bottom, the Her (V = 186.9 kms"'), the Gau {V = 183.2 kms"'), 
theMoml (V = 186.2 kms"'), the Bulk (V = 186.1 kms"') and the T(peak) velocity fields (V = 185.4 kms"'). The curve shows the Hermite fit. Right: Observed 
profile of the Hi emission in IC 2574 at position (o, (5)(2000.0) = 10*28"'51.4',68°28'13.0" (thick line). The five vertical lines indicate the derived typical 
velocities using, from top to bottom, the Her (V = 64.2 kms"'), the Gau (V = 66.3.2 kms"'), the Moml (V = 68.1.2 kms"'), the Bulk (V = 80.8 kms"') and the 
T(peak) velocity field (V = 63.8 kms"'). 



An additional source of uncertainty is formed by possible 
differences in r otation velocity between the approaching and 
receding sides. ISwatersI (Il999h assumes that the differences 
between the rotation curves derived for the entire disk and 
those for either the approaching or receding side represent 
a 2(7 difference, and defines a pseudo-lcr uncertainty due to 
asymmetries, as one fourth of the difference between the ap- 
proaching and receding side velocities. Though this is only an 
assumption, we follow this convention and define the uncer- 
tainties in the rotation curves as the quadratic addition of the 
dispersions found along the rings and the pseudo-la uncer- 
tainties due to asymmetries between approaching and reced- 
ing sides. Note that this is a very conservative definition: the 
difference in velocity between approaching and receding sides 
will already be partly reflected in the dispersion found along 
a ring. The uncertainties in the rotation velocities presented 
here thus give a realistic, and possibly even conservative, pic- 
ture of the error budget of our rotation curve determinations. 

4. INDIVIDUAL GALAXIES 

Detailed descriptions of the data as well as the derivation 
of the tilted ring models and the resulting rotation curves of 
individual galaxies are given in the Appendix. In this section 
we compare these results with previous determinations of the 
various Hi rotation curves from the literature. (For NGC 362 1 
and NGC 3627 we could find no previous, published determi- 
nations of their H I rotation curves. These are therefore not 
discussed in this Section and we refer to the Appendix for 
more information.) We do not compare with Ha or CO ro- 
tation curves. These are determined using different methods, 
and originate from different phases of the ISM; systematic ef- 
fects due to these differences are not straight-forward to quan- 
tify and beyond the scope of the current paper 

For reference, a summary of the resulting tilted-ring pa- 
rameters is given in Table |2] These include the position 



of the dynamical center, the rotation curve sampling inter- 
val, systemic velocity and radially averaged values of PA 
an d i. For the central positi ons we use the results presented 
in iTrachternach et all (l2008h as derived using th e same data. 
We ad opt luminosities and distances as given in I Walter et alJ 
(1 2008 !). and where necessary literature measurements have 
been corrected to the same distance scale. 

4.1. NGC 925 

NGC 925 is classified as a late-type bar r ed spi ral. It was 
previously observed in H I by iPisano et al.l ( Il998h . Figure |6] 
compares the rotation curves. The main difference with our 
analysis, namely the position of the dynamical center, is dis- 
cussed in ITrachternach et"aLl ( l2008l) . This difference also ex- 
plains the slightly higher value for Vsys of 551.5 kms"' found 
bv lPisano et al.l ( Il998h . compared to our 546.3 kms"'. Other 
than this, we derive broadly similar trends of ; and PA with 
radius. The curves themselves agree within the uncertainties, 
showing that for solid-body (inner) rotation curves, such as 
this one, the shape does not critically de pend on the position 
of the center. The uncertainties quoted in lPisano et al.l (1 19981) 
are the formal fit errors. As discussed before, these are un- 
realistically small, representing the error in the mean velocity 
along a ring, rather than the scatter In this case these formal 
errors are smaller than the sizes of the symbols in Fig.|6] The 
uncertainties we have adopted are more representative of the 
true physical uncertainties. 

4.2. NGC 2366 

The data and analysis of NGC 2366 are given in lOh et alJ 
(l2008h . We refer to their paper for a complete discussion. 

4.3. NGC 2403 
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Fig. 4. — Rotation curves forNGC 2403 derived using the five different velocity fields. Dotted curve: Moml; dashed curve: Gau; full curve: Her; long-dashed 
curve: Bulk; dot-short dash: T(peak). Also shown are the differences with respect to the Hermite curve. Thick black curve: Her-Gau; thick gray full curve: 
Her-Moml; thick dotted gray curve: Her-Bulk; thick dashed gray curve: Her-T(peak). The differences have been multiplied by a factor of 5 for clarity. Bottom: 
Rotation curves and differences for IC2574. Line types as described above. 



NGC 2403 is a late-type Sc spiral and member of the M8 1 
group. Its H I rotation curve has been derived in many pre- 
vious studies. The firs t measur e ment using a synthesis tele- 
scope was presented in lShostafl (Il973h . NGC 2403 has since 
been re-observed and reanalyzed at ever increasing sensitivity 
and re s olution by amo ngst others Bosma ('1978'), iBegemanI 
(Il987h . ISickind Jl997h and iF^aternali et al (2002i). In Fig. |7] 
we compare our curve with some of these analyses, all cor- 
rected to the same distance. There is generally good agree- 
ment between the various determinations. Beyond ~ 13 kpc 
the curves differ somewhat, but most of the velocity infor- 
mation at these radii originates from close to the minor axis 
of the velocity field, and thus depends on details such as the 
weighting used in the tilted ring fit, the flux limit used to mask 



the velocity field, etc. 

Other parameters also compare fav orably: the average / and 
PA values fo und bvlBegemanI d 19871) were 60.2° and 122.5°, 
respectively. ISicking|(ll997h found 60.9° and 123.9°, whereas 
iFraternah et all (l2002l) derived 62.9° and 124.5°. Our aver- 
age values (defined as the average unweighted values of the 
model distributions) agree well, at 62.9° and 123.7°, respec- 
tively. Similarly, the average absolute difference between our 
Vsys value and the three literature values is only 0.5 kms~'. 

The "bump" at 7? ~ 15 kpc in the THINGS curve is real, in 
the sense that it cannot be explained by an incorrectly chosen 
/ or PA. Most of the information in this part of the curve orig- 
inates from the two arms that are seen to stretch away in the 
outermost northern and southern parts of the H I map. Mo- 
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Fig. 5. — Top panel: Comparison of input model rotation curves with "observed" curves as derived from model data cubes smoothed to the THINGS 
resolution. The full curve represents the input model rotation curve of a model that rises to its maximum rotation velocity within 1 kpc. The open circles 
represent the corresponding rotation curve as observed at the THINGS resolution. The dashed curve and open squares are the corresponding curves for a model 
rising to its maximum value within 0.5 kpc. Bottom panel: Difference between "observed" rotation velocity and input model velocity. The full line represents 
the 1 kpc model; the dashed curve the 0.5 kpc model. 



tions along these arms could explain the feature, as suggested 
by the kinks in the velocity contours at those radii. The veloc- 
ity fields of NGC 2403 are also discussed in Sec. 3.2. 

4.4. NGC 2841 

NGC 2841 is an early-type (Sb) spiral and was previously 
observed at lower resolution in H I by Begeman ( 1987). We 
compar e the curves in F ig. [8] and find good correspondence 
overall. iBegemaril ( Il987h finds a systemic velocity of 63 1 . 1 ± 
1.4 kms"', a difference of only 2.6 kms"' with our value of 
633.7 lb 1 . 8 krn s"'. The average PA and / values found by 
lBegemanl ( ll987h are 154.5° and 73. 1°, respectively, and differ 
only by a few degrees from our values (cf. Table|2li. 



4.5. NGC 2903 

NGC 2903 was previously observed in H I by iBegemanI 
(119871). Fig. |9] compares the two rotation curves. The 
IBegemanI (HgST) value of V,y, = 557.3 ± 1.3 kms"' agrees 
well with our value of 555.6 ± 1.3 kms"'. The outer de- 
clining parts of the rotation curves also agree well with each 
other There are however significant differences in the inner, 
rising parts of the rotation curve, where the curve by lBegemanl 
11987,) rises much more steeply. His data did not allow a full 
tilted-ring model in the inner part of the galaxy, and the sud- 
den changes in PA and / are therefore also not dete cted in his 
analys is. Rather, the innermost data points of the IBegemanI 
( Il987l) curve were derived from a major axis position-velocity 
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Fig. 6. — Compa rison of the NGC 92 5 THINGS rotation curve with the curve derived in IPisano et alj <1998r) . Symbols as indicated in the figure. The 
uncertainties given in IPisano et alj jl998l) are smaller than the symbol sizes. However, there uncertainties are the formal least-square fit errors and severely 
underestimate the true uncertainties in the measured rotation velocity. Cf. discussion in Sect. 3.5. 



diagram. The higher spatial and velocity resolution of our 
data and the different methods used to extract the velocities 
both contribute to this large difference. 

4.6. NGC 2976 

NGC 2976 is cl assified as an Sc gal axy and was previously 
observed in H 1 bv lStil & Israeli ( |2002|) . They chose the center 
of mass of the H I disk as their dynamical center, and there is 
an offset of ^ 10" with respect to our choice. Their choice for 
Vsys is also slightly different from ours: 4±2 kms"', versus 
ourl.l±1.3 km s-i . |StiL&Israel| (gOOl) find ; = 65 ° and PA 
= 326°, versus our values of 64.5° and 334.5°, respectively . 
Figure [To] compares the two curves. The IStil& Israeli (I2OO2I) 
curve suggests a flattening towards the outer radii, whereas 
our curve keeps rising until the last measured point. The 
curves do, however, agree within their respective uncertain- 
ties. 



4.7. NGC 3031 

Better known as M8 1 , this galaxy is the proto-typical grand- 
design spiral and togeth er with M82 and NGC3077 forms a n 
interacting system, (e.g.. Ivan der HulsjflQTliYun et a]lll994ll . 
A previous determ ination of the Hi rotation curve was by 
IVisseil(ll978lll980l) . We compare our results in Fig. [TT] after 
correction to the same distance scale. We find excellent agree- 
ment overall between his and our rotation curves where they 
overlap in radius. The differences we find are all attributable 
to different local values for PA and we allow ; and PA to 
change with radius, whereas fVisserl (Il980l) uses a constant in- 
clination of 59° and a PA value of 329°. These constant val- 
ues agree, however, well with our average values of 59° and 
330.2°, respectively. 



4.8. NGC 3198 
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Fig. 7. — Comparison of the NGC 2403 THINGS rotation cui-ve with recent determinations from the literature. Symbols and references are indicated in the 
figure. 



NGC 3 198 is known for its prot o-typical flat rotation curve 
(see e.g. Ivan Albada et al.l [19851) . We compare our rota- 
tion curve with tho se derived by iBegemaril (119871 119891) 
and ISickind ( 119971) in Fig. [121 The systemic velocities 
agree within their res pective loi uncertainties: our result is 
660.7±2.6 kms "'.rBegemani (ri98a) finds 660.4±0.8 kms"', 
whereas [Sickingl (Il997l) finds 659.4 ±2.6 km s The respec- 
tive average PA and ; values only show differences < 1 ° . 

Turning now to the rotation curves themselves and starting 
with the outer parts, we find that our curve agrees well with 
the Sick ing ( 1997) result, but is a few kms"' higher than the 
iBegemanI ([1989) data. Noting that the uncertainties in neither 
the Sicking nor the Begeman analysis do take into account 



the difference between approaching and receding sides, the 
results are still consisten t with each other. 

In the inner parts, the IBegemaril (119891) data show consis- 
tently higher rotation velocities. These were, however, de- 
rived from major-axis position velocity diagrams, rather than 
from a tilted-ring model. Due to the much lower resolution 
they were also explicitly corrected for beam smearing. The 
sys tematic offset co uld therefore imply an overcorrection of 
the lBegemanl ( Il989l) data, or a system atic deviat i on due to the 
use of major axis profile velocities. IBegemaril ( 119871) men- 
tions that the beam smearing corrections are only significant 
within 1.5' 6 kpc), precisely the radius within which the 
curves show disagreement. The inner parts of our curve are 
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Fig. 8. — Compaiison of the NGC 2841 THINGS rotati on cur ve with the analysis by Begeman ( 1987). Symbols and references are indicated in the figure. 
Note the different nature of the uncertainties: the Beg ernMil j 19871) result uses the formal errors in the least-squares fit, whereas our definition takes into account 
the differences between the two sides of the galaxy, and the dispersion in velocity values found along the rings. 

fully consistent with the lSickin^(ll997h rotation curve. 
4.9. IC2574 

The data and analysis of IC 2574 are described by lOh et aP 



(l2008|)and we refer to their paper for a more complete discus 
sion. The velocity fields of IC 2574 have also been discussed 
in Sec. 3.3. 

4.10. NGC 3521 

NGC 3521 was pr e viousl y observed in Hi by 
ICasertano & van GorkomI (Il991h . Its rotation curve was 
presented as one of the first cases where a genuinely de- 
clining H I rotation curve was found in a spiral galaxy. We 



compare our curve with the ICasertano & van GorkomI (1199 Ih 
curve in Fig. [T3] The agreement is not perfect, the literature 
curve underestimates the rotation velocity in the inner parts, 
whereas our curve d oes not show the steep drop in th e outer 
parts. Unfortunately. ICasertano & van GorkomI (1 19911) do not 
list the values they assume for / and PA, which precludes 
a more detailed comparison. The high inclinations, and 
resulting insensitivity of the rotation velocities to the exact 
inclination values, make it unlikely that the difference is 
a pure inclination issue. For example, a change in ; of 5 
degrees, starting from an initial value of / = 75°, only amounts 
to a change of some 2 percent or ~ 4 kms~', not the ~ 30 
kms"' that is needed to bring the ICasertano & van GorkomI 
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Fig. 9. — Comparison of the NGC 2903 THINGS rotation curve witli tlie analysis bv lBegemanl 119871) . Symbols and references are indicated in the figure. See 
Fig.[H]for a description of the uncertainties. 



(1199 Ih in agreement with ours. In summary, even though we 
find rotation velocities in the outer parts that are lower than 
in the inner parts, we do not find evidence for the steep drop, 
n or for the negative gradien t in ou ter velocities as presented 
in lCasertano & van Gorkoml(ll991h . 

4.11. NGC 4736 

NGC 4736 was previously observed in Hi by iMuldej 
(Il995h who also derived a rotation curve. Fig. [14] shows a 
comparison between the respective rotation curves. Due to 
the lower resolution and sensitivity of his data, [Mulder ( 1995|) 
did not detect a trend of inclination with radius, and therefore 
assumed a constant value of 40°. In Fig. [14] we therefore also 



show his rotation curve corrected using our incUnation val- 
ues. Although the large-scale trends are similar to those found 
for our curve, the THINGS curve in general shows more pro- 
nounced small-scale features. 



4.12. DD0154 

DDO 154 is a gas-rich dwarf galaxy whose rotation curve 
has been extensively stud ied (see lCarignan & Freeman|[T988l : 
ICarignan & Purton]|1998h . It was also was one of the first 
galaxies whose rotation curve was used to illustrate t he con- 
flict with Cold Dark Matter predictions (iMoordl 1994b . DDO 
154 shows clear evidence for a warp in the outer parts. There 
is no evidence for large non-circular motions in the veloc- 
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Fig. 10. — Comparison of the NGC 2976 T HINGS rotation curv e with the analysis by Stil & Israel (200?). Symbols and references are indicated in the figure. 
Note the different nature of the uncertainties: IStil & Israel 120021) use the difference between approaching and receding side as a measure for the uncertainty, 
whereas our definition additionally takes into account the dispersion in velocity values found along the rings. 

Given the differences found in the tilted ring parameters 
of the approaching and receding sides separately , the er- 
ror bars belonging to the ICarignan & PurtonI (Il998h rotation 
curve most likely underestimate the true uncertainties. Also 
note that our rotation curve in general shows h igher rota- 
tion velocities than the ICarignan & f^irtonl (Il998h one. This 
is not merely an inclination effect, as our values for / are 
in general higher by a few d egrees than the ones presented 
in ICarignan & f^irtoiil (Il998h . Correcting the curves to the 
same inclinations would therefore only increase the differ- 
ence. Note that a model with a decreasing inclination in the 
outer parts, as perhaps suggested by our analysis (cf. the Ap- 
pendix), would lead to an increasing rotation velocity in the 



ity field (iTrachternach et al"1l2008h . This is evidenced by the 
residual velocity field where the absolute value of the residu- 
als nowhere exceeds ~ 6 kms"'. 

Previous d eterminations of the rotation curve of DDO 
154 ( see e.g. . ICarignan & Freemanl[l988HCarignan & PurtonI 
Il998h . show an apparent drop in the rotation velocity in the 
outermost parts, which was interpreted as evidence that at 
these radii one has almost reached the edge of the dark mat- 
ter halo. In Fig. [15] we compare our rota tion curve with one 
of the se det erminations. Contrary to the ICarignan & PurtonI 
(119981) and ICarignan & Freeman! (Il988l) results, we do not 
find strong evidence for a decline and our rotation curve in 
the outer parts is consistent with being flat. 
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outer parts. The only way in which our determination of the 
rotation curve could be made consistent with a declining rota- 
tion curve is to increase the inclination, which is ruled out by 
our data. 



4.13. NGC 4826 

NGC 4826, or the "Evil-Eye" galaxy is an early-type spiral 
with as its most remarkable kiner natical feature the presence 
of two counter-rotating gas disks jBraun et al.lll994l) : the in- 
ner high-column density H 1 disk is associated with the bright 
stellar disk and rotates in the same direction. The outer, much 
lower column-density H I disk rotates in the direction opposite 
to that of the stars. The transition between the two disks oc- 



curs around ^ 100" radius. iBraun et"an ( 119941) speculate that 
the outer disk might find its origin in a merger of a gas-poor 
spiral with a gas-rich and star-poor dwarf with a retrograde 
spin (with the mass ratio of the dwarf and the spiral < 0.1; 
iRix et al.lfT995h . 

In Fig. [T6l we c ompare our rotation curve with that from 
iBraun et al.l Jl994l) . Even though to first order there is rea- 
sonable agreement, the scatter in the velocities is much larger 
than in any of the other rotation curves presented here. Braun 
et al. remark that between 100" <R< 300" (4 kpc < R < 
12 kpc) the kinematical parameters are difficult to determine 
due to what they call a "severe distortion" in the velocity field, 
showing that the transition region between the two counter- 
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Fig. 12. — Left: Comparison of the NGC 3198 THINGS rotation curve with recent determinations from the literature. All curves have been corrected to a 
distance of 13.80 Mpc. Symbols and references are indicated in the figure. 



rotating disks deviates from circular moti ons. Similar conc lu- 
sions were reached for the ionized gas bv lRix et"aLl d 19951) . 

4.14. NGC 5055 

NGC 5055 is an Sbc galaxy with an extended and warped 
tenuous outer H I disk. It s Hi distributio n and kinematics 
were recently analyzed by iBattaglia et al.l (12006,) . For Vsys, 
iBattaglia et all (|2006|) find a value of 497.6 ±4.8 kms"', 
which is in close agreement with our value of 496.8 ±0.7 
kms~'. Their distributions of / and PA are also in agreement, 
and they note a similar difference between the models for ap- 
proaching and receding sides. 

We compare the two rotation curves in Fig. [17] The first 
inner maximum at /? ~ 2 kpc is more pronounced in our data 
set. This is likely due to a combination of increased resolution 
and the use of Her mite polynomia l s as op posed to the Gaus- 
sian functions that IBattaglia et al.l (|2006|) used. This is also 
the most likely explanation for the higher rotation velocities 



we find in the outer parts. 

4.15. NGC 6946 

NGC 6946 is a late-type spiral which, in terms of H I stud- 
ies, is better known for its population of H I structures such 
as H i holes and high-v elocity gas (e.g.. lKamphuis & Sancisil 
119931: lBoomsmall2007h . Its global dynamics have been rela- 
tively little studied. This can be partly attributed to its fairly 
low inclination which is at the limit of what is feasible using 
tilted-ring studies. 

The rotation curve of NGC 6946 was determi ned before 
b y ICarignan et al.l (Il990h and more recently by iBoomsmal 
(200 7|. The comparison w ith our curve is made in Fig. [Tsl 
Both lCarignan et al.l (1 19901) and lBoomsinal (l2007h used a con- 
stant value of /= 38° in their analyses and these different incli- 
nations, compared to our average value of 32.6°, explain the 
offset between the curves. Apart from thi s, the only promi- 
nent difference is the large rotation velocitv lBoomsmal ( l2007h 
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Fig. 13. — Compai'ison of the NGC 3521 THINGS rotation curve with the analysis bv lCasertano & van Gorkon]d 19911) . Symbols and references ai'e indicated 
in the figure. 



finds for the innermost point of the rotation curve. In our data 
the corresponding area of the velocity field is below our 3a 
H I column density level, and we can therefore not comment 
on the presence of this feature in our data, though there is a 
hint of a turn-up in the major axis position-velocity diagram 
(cf. Fig. 



ity resolutions and different choices for the inclination distri- 
bution (see dis cussion in the Ap pendix). Note that the inner 
8 points of the lBegemanI (119871) curve are based on position- 
velocity diagrams and not on a tilted-ring fit to a velocity field. 
The differences in the outermost parts are due to the different 
inclinations used there. 



4.16. NGC 7331 

NGC 7331 is an early-type Sb spiral with prominent spiral 
arms. It was previously olDserved in H 1 by [Begeman (1987]). 
Figure[T9]compares both curves. It is clear that the two curves 
broadly agree, but there are many small differences. Most of 
them can be attributed to the very different spatial and veloc- 



4.17. NGC 7793 

NGC 7793 is a late-type Sd spiral gal axy in the Sculptor 
roup and was observed before in H 1 by ICarignan & Puchg 
1990). In Fig.|20]we compare our rotation curve with their 
results. The dyna mical center and systemic velocity, as well 
as PA, as found bv lCarignan & Puchd ( Il990h agree with ours 
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Fig. 14. — Com parison of t he NGC 4736 THINGS rotation curve (grey points) witli the analysis by Mulder (1995). The open circles show the rotation curve 
as published in .Muldej j 19951) assuming a constant inclination of 40° . The filled black circles show the same data but corrected using our inclinations. 



within the uncertainties. Th e behavior of the inclinat ion is 
subtly different: whereas the lCarignan & Puchd (Il990h anal- 
ysis indicates a rising inclination in the outer parts, we find a 
(better defined) gradual decrease there. The resulting rotation 
curves on the whole agree reasonably well with each other, 
but differ in small but important respects. The most important 
one is the more gradual decline of our rotation cu rve in the 
outer parts. Even though ICarignan & f^ichd (Il990h note that 
this is "one of the very few cases with a truly declining rota- 
tion curve" (their emphasis), we argue that any drop in rota- 
tion velocity is much less extreme, especially as there is some 
evidence that the inclination might decrease more strongly in 
the outer parts than we have assumed (cf. the distribution of 



inclination with radius of the approaching and receding sides 
in Fig. [86b . Adopting this steeper drop would raise the ro- 
tation velocities we find in the outer parts by ^ 25 kms"'. 
We stress that these inclination values are uncertain, and de- 
termined by only a small number of pixels, but they seem to 
suggest that a declining rotation curve cannot be unambigu- 
ously established. 

4.18. Summary 

In summary, we find in general good agreement with pre- 
vious determinations of the rotation curves of the galaxies in 
our sample. The most notable differences are: 

• We find some differences with previous determinations 
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Fig. 15. — Comparison of the DDO 154 THINGS rotation curve with the curve derived in lCarignan & PurtonI 1199811 . Symbols and references are indicated in 
the figure. 



of the inner rotation curves of NGC 2903 and NGC 
3198. Both galaxies contain a small central bar. The 
presence of this feature combined with the lower reso- 
lution of the literature data can contribute to this differ- 
ence. 

• The outer rotation curve of NGC 5055 is higher than 
previously determined. Most likely our use of a Her- 
mite velocity field contributes to this difference. 

• We do not find steep declines in velocity in the outer 
rotation curves of NGC 3521, NG C 7793, DDO 154 
and NGC 2366 (see lOh et al1l2008l for the latter curve). 
Where declines are observed, they are more gentle, and 
(within the uncertainties in rotation velocity and incli- 
nation) consistent with flat rotation curves. 



5. MASS MODELS 



We now use the THINGS rotation curves, in combination 
with information on the distribution of gas and stars, to con- 
struct mass models of our sample galaxies. These models are 
then used to quantify the distribution of dark matter within 
the galaxies. We calculate our best model estimate for the 
observed velocity as 



^ obs "^eas^^*"^*^"^! 



2 

halo" 



(3) 



The mass models thus need as input the observed rotation 
curve, Vobs, the rotation curve of the gas component, Vgas, and 
the rotation curve of the stellar component, y*. The rotation 
curve from the halo, Vhaio, can then be calculated from these 
known input curves. An additional free parameter is the stellar 
mass-to-light ratio, T*, which is introduced as we generally 
can only measure the distribution of the light of the stellar 
population, rather than the required mass. This parameter in- 
troduces one of the largest uncertainties in the procedure and 
will be discussed in detail in Sect. 5.2. 
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Fig. 16. — Compai'ison of the NGC 4826 THINGS rotation curve with the analysis by Braun et^ <1994ft . Symbols as indicated in the figure. Plotted are the 
absolute rotation velocities, the inner disk (i? < 4 kpc; indicated by the horizontal dotted line) rotates in the opposite direction compared to the outer disk. 



While the concept of mass modeling is well described by 
Eq. [3] there are some intricacies that warrant a brief discus- 
sion. The rotation velocities under consideration are in all 
cases the velocities of respective test particles moving in cir- 
cular orbits in the plane of the galaxy. The rotation velocities 
of the gas, Vgas, and the stars, V^,, are defined as the velocities 
that each of the components would induce on a test particle 
in the plane of the galaxy, if these disks were sitting in isola- 
tion, without any external influences. These velocities in the 
plane are calculated from the observed baryonic mass density 
distributions, taking into account the (assumed) vertical den- 
sity distribution. For the gas disk one generally assumes a 
thin disk, while for the stellar disk a thick, sech" distribution 



is chosen. These procedures and choices of parameters are 
described in more detail on Sect. 5.1 for the rotation curve of 
the gas and in Sect. 5.2 for that of the stars. The latter section 
also describes in detail possible choices for T*. 

The halo velocities, Vhaio, in Eq. [3] (or rather the velocities 
induced on a test particle by the halo potential) are also mea- 
sured in the plane of the galaxy and assume a spherical halo. 
Sect. 5.3 briefly discusses the dark matter models used, and 
in Sect. 5.4 the procedure of constructing the models is pre- 
sented. Specific discussions dealing with individual galaxies 
are given in Sect. 6. 

Initial analysis of the rotation curves of NGC 3627 and 
NGC 4826 showed their dynamics to be severely affected by 
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Fig. 17. — Compaiison of the NGC 5055 THINGS rotation curve with the analysis bv lBattaglia et al.l j2006l) . Symbols and references are indicated in the 
figure. 



non-circular motions (iTrachternach et al"]l2008h and we will 
therefore not consider these two galaxies further. 

5.1. Neutral Gas Distribution 

To derive the rotation curve associated with the neutral gas 
component, we use th e primary beam and residual-scaled in- 
tegrated Hi maps from lWalter et alj (120081) . We derive neutral 
hydrogen surface density profiles using our previously derived 
tilted ring geometrical parameters. The surface densities are 
corrected by a factor 1 .4 to take into account the presence of 
helium and metals (associated with the atomic hydrogen; the 
helium component by itself scales with a factor of 1.33). The 
ROTMOD task in GIPSY is then used to compute the corre- 



sponding rotation curve, under the assumption of an infinitely 
thin gas disk. 

As we are dealing with a mass distribution in the form of a 
disk, rather than a sphere, the rules-of-thumb that apply to the 
gravitational force of a spherical mass distribution do not ap- 
ply. Most importantly, the convenient result that the mass dis- 
tribution outside a particular radius does not affect the effec- 
tive force inside that radius is not universally valid for a disk. 
Specifically, in the case of a disk with a sufficiently promi- 
nent central depression in the mass distribution, a test particle 
located in or near it can feel a net outward force. A literal 
interpretation in the context of normal gravity would imply 
negative mass, an imaginary rotation velocity and therefore a 
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Fig. 18. — Comparison of the NGC 6946 THINGS rotation curve with the analyses bv lCarignan et al.l I199CD and lBoomsmal 120071) . Symbols and references 
are indicated in the figure. 



negative value of Here, following the convention in rota- 
tion curve literature, we plot the imaginary rotation velocities 
as negative velocities, keeping in mind that it really is V~ that 
is negative. We stress that this convention for representing the 
imaginary velocities is merely adopted for convenience and 
does not imply counter rotation, nor the presence of repulsive 
matter or other exotic phenomena: it simply is a short-hand to 
enable illustrating the effective force on a test particle caused 
by a non-spherical mass distribution with a central minimum. 
The mass model of NGC 3031 in Sect. 6.6 is a good example 
of this phenomenon. 

5.2. Stellar Distribution 



One of the main uncertainties in deriving mass models is the 
conversion from observed stellar luminosity to stellar mass 
(through the stellar mass-to-light ratio T*). In the optical 
bands, is dominated by the young stellar population, recent 
star formation events and the effects of dust and/or metallic- 
ity. Values of derived in, e.g., the B-band thus have in- 
herently a large uncertainty which can be up to a factor of 
a few. These limitations are the reason why extreme mass- 
models involving so-called " maximum" or "minimum disks" 
(Ivan Albada & S and sil [19861) are often used. The maximum 
disk assumption maximizes the contribution of the baryonic 
matter by assuming the maximum value of allowed by the 
total rotation curve. Minimum disk is the opposite and mini- 
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Fig. 19. — Comparison of the NGC 733 1 THINGS rotatio n curve with result bv lBegemanI 1198711 . Symbols and references as indicated in the figure. Note 
the different nature of the uncertainties: the lBegemanI jl987l) result uses the formal errors in the least-squares fit, whereas our definition talces into account the 
differences between the two sides of the galaxy, and the dispersion in velocity values found along the rings. 

mizes the baryonic contribution by setting T,^ = 0. Sometimes 
the contribution of the gas is also ignored (and all rotation 
is therefore assumed to be induced by the dark matter halo). 
These two extreme T,^ assumptions can thus be regarded as 
"minimum halo" and "maximum halo" models, respectively, 
minimizing or maximizing the contribution of the dark matter 
Though not necessarily physically plausible, they do provide 
hard upper and lower limits on the concentration and distribu- 
tion of the dark matter 

Fortunately, the uncertainties in decrease dramatically 
towards the near-infr ared (from a factor ^ 10 in the B-band to 
a factor ^3mK; cf. lBell&de Jongll2001h . At near-IR wave- 
lengths one mostly probes the old stellar pop ulation contain- 
ing the bulk of the stellar mass, as described in lBell & de Jond 
(|200lj) and Verheijen (1997). The range in plausible values 
for is in the near-IR only a shallow function of the color of 
the disk, and observations at these wavelengths are therefore 



useful in constraining the mass of the disk. IVerheijenI (Il997h 
finds that maximum disk fits to a sample of HSB galaxies in 
the UMa c luster require only a srn all range of Tf =0.7 ±0.2. 
Similarly, iBell & de Jong! (1200 1 h find, from population syn- 
thesis modeling, that disks of galaxies exhibit a typical range 
of Tf = 0.5 ± 0.2 measured over the full range of observed 
galaxy colors. 

To characterize the stellar distributions in o ur sample galax- 
ies, w e use the 3.6 fim images from SINGS (iKennicutt et al.l 
l2003h . (For NGC 2366, NGC 2903 and DDO 154, which are 
not part of the SINGS sample, we have retrieved data from the 
Spitzer archive.) These pro vide a good proxy for the emission 
of the (old) stellar disk (cf. lPahre et al.ll20o4 . Note that some 
of the 3.6 fim emission may originate in AGB stars, PAHs and 
hot young stars, but as we will only be using the radial distri- 
butions of the emission, their contribution (if present) will be 
strongly diminished by the azimuthal averaging. 
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TABLE 2 

TILTED-RING MODEL PARAMETERS OF THE THINGS GALAXIES 
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3.9 
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226.2 


49.6 


290.1 



Note. — (1) Name of galaxy; (2) Right Ascension (J2000.0); (3) Declination 
(J2000.0), center p ositions from Trachtemach et al. (2008); (4) Distance as given in 
[Walter et alj 120081) ; (5) Sampling increment used to derive the rotation curve; (6) 
Adopted systemic velocity; (7) Average value of the inclination; (8) Average value of 
the position angle of the receding side, measured from north to east, and in the plane 
of the sky. 



After the images were cleaned of foreground stars, the tilted 
ring parameters were used to derive the radial luminosity pro- 
files. We used the standard photometric calibration provided 
with the SINGS images to construct the surface brightness 
profiles. In most cases the 3.6 /im surface brightness distri- 
butions can be well described by a single exponential disk. 
In a small number of galaxies, an additional central compo- 
nent was seen, containing only a small fraction of the total 
luminosity of the galaxy. It is a matter of debate whether it 
is appropriate to call this central component a "bulge". The 
surface brightness distributions measured are not as steep as 
the canonical bulge surface distribution, which may re- 
flect the THINGS sample selection criteria, with its emphasis 
on more late-type galaxies. Beyond this, the radial range over 
which this inner profile dominates the total emission is not 
large enough to determine accurately enough the value of n in 
an parameterization of the surface brightness profile of 
that central component. This is, however, not a major prob- 
lem due to its rapid drop in surface brightness and we found 
that most often an exponential surface brightness distribution 
provides an adequate and convenient fit. We experimented 
with a number of different functional power-law forms, but 
found that this did not affect the shape of the (compact) rota- 
tion curve of this component to any relevant degree, and does 
not impact on our final mass models. Where necessary, we 
therefore decomposed the surface brightness profile in inner 
and outer exponential components. In the rest of this anal- 
ysis we refer to the outer component as the "outer disk", or 
simply "disk". The inner component we interchangeably call 
"inner disk" or "bulge". This simply is a short-hand to indi- 
cate the components, and, especially for the inner component, 
should not be taken as a statement on its implied evolution or 
other physical properties. Central surface brightness values 
and scale lengths, where mentioned in the text, are not cor- 
rected for inclination. Needless to say that this correction has 



been applied prior to calculation of the rotation curve. 

As far as we are aware, no exhaustive investigation has yet 
been made into the relation between 3.6 /im emission and T,^. 
We therefore use an empirical approach based on results de- 
rived using the 2MASS J, H and K bands to derive T* in 
the 3.6 /im band (hereafter T^-^) from the Spitzer observa- 
tions. All our sample galaxies are included in the 2MASS 
Large Galaxy Atlas (^larrett et al 2003i), and we use the J-K 
colors given ther e, together with the re lations between color 
and Ti, given in iBell & de Jond (1200 Ih. as we ll as the rela- 
tion between T^^ and T^^ from loh et all (l2008l) . to derive the 
relation between T^^ and J-K. We find 

log(Tf)= 1.43(7-/*:)- 1.38 (4) 

and 

Tj'' = 0.92 Tf- 0.05. (5) 

One of the largest uncertainties in determining T^, is the 
stellar Initial Mass Function (IMF). Some asp ects of this, 
within the cont ext of THINGS, are di scuss ed in [ Lerov et al] 
(I2008h : see also lBell&de Jong! (1200 Ih and lBell et al.l J2003h 
for a more gen eral discussion. Using stellar popula tion syn- 
thesis models, iBell & de Jond ( 1200 II) show that a ISalpeteij 
(1955) IMF results in stellar disk masses that are too heavy 
to be consistent with dynamical (maximum disk) constraints 
from rotation curves. By scaling the stellar disk masses down 
by a factor of 0.7, in order to make them consistent with the 
rotation curves, they introduce a "dief'-Salpeter IMF, equiva- 
lent to a standard Salpeter IMF, but with a reduced number of 
stars below 0.35 Mq. This diet-Salpeter is thus the maximum 
IMF (i.e., yielding the highest disk masses for a given photo- 
metric constraint) allowed by galaxy dynamics studies. Stud- 
ies of the stellar population in the Milky Way suggest an IMF 
that produces lower disk masses (iKroup a 2001). The diet- 
Salpeter disk masses are 0.15 dex (a factor 1 .4) rn ore massive 
than the Kroupa ones. We refer to iLerov et aTl (|2008.) and 
iBell et all ( l2003h for a discussion on other IMF indicators. 
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The disk masses as impUed by Eq. (4) assume the diet- 
Salpeter IMF, and are therefore the maximum disk masses 
allowed for a given J -K color. To gauge the effects of a 
"lighter" IMF on our rotation curve results we also derive disk 
masses assuming a Kroupa IMF by decreasing the constant 
term in Eq. (4) by 0.15 dex. 

Disk galaxies typically show radial color gradients of a few 
tenths of a magnitude, becoming bluer with radius. This is 
interpreted as a combined stellar age and metallicity gradient, 
with the outer regions younger and of lower metallicity; dus t 
seems to play only a minor role in this regard dde Jondll996l) . 
In the optical, these gradients are thought to imply radial 
differences of a factor of ~ 2 between inner and outer parts 
of galaxies. Even in the /iT-band where many of the effects of 
recent star formation and dust are minimized, changes in T^, 
of ~ 20 to 30 per cent are not uncommon. We therefore 
associate the (J-K) color gradients observe d in the disks of 
our sample galaxies with variations in T*. lOh et al.l ( l2008h 
show that similar gradients in T^^ are expected in the 3.6 fim 
band'^. We refer to Ide JongI (1996) for an extensive discus- 
sion of color gradients in the optical and near-IR. Our goal 
is to create mass models taking into account radial varia- 
tions within the stellar disk. We determine "T^^-weighted" 
mass-surface density profiles for our galaxies and use these as 
input for ROTMOD to compute the rotation curve of the stellar 
disk. For these stellar disks we assume a vertical sech^ distri- 
bution and use zo = h/5 where h is the exponential radial scale 
length ( Ivan der Kruit & SearlellI98Tallbl) '^. The rotation curve 
derived in this way does not critically depend on the value of 
h/zo- 

One component we have not explicitly addressed is the 
molecular gas. However, as the molecular gas is (at least in 
non-dwarf) THIN GS disk galaxies distributed like the stars 
(iLeroy et al.ll2008h . and with a surface den sity only a few per- 
cent of that of the stars (iPortas et al.ll2008l) . its (small) contri- 
bution is reflected in a small i ncrease in T< ,. 

In an accompanying paper, ILeroy et al.l (l2008h present an 
alternative empirical method to derive values from the 3.6 
images. They use the 3.6 fim and 2MASS K-hand im- 
ages themselves to directly derive an empirical pixel-to-pixel 
relation between the two bands, enabling them to predict K- 
band fluxes from the IRAC images. Combining these /iT-band 
predictions with a constant Tf = 0.5 then yields the stellar 
masses. We compare with our results (discussed in Sect. 6) 
in Fig. [21] for both our assumptions for the IMF. The masses 
we derive tend to be slightly larger, a least-squares fit yields a 
slope of 0.94 with a scatter of 0.12 dex. Given the empirical 
nature of both methods, and the absence of a full set of rigor- 
ously computed population synthesis models, this should be 
regarded as satisfactory agreement. 

5.3. Dark matter halos 

Here we describe the two well-known models for the dark 
matter distribu tion that we will use in our analysis . Firstly, 
the NEW halo (iNavarro. Frenk & WhitdfT996lfT997h . derived 
from CDM simulations, and dominated by a central density 
cusp. Secondly, we discuss the venerable, observationally 

Note that, as mentioned before, a small amount of dust emission due 
to, e.g., PAHs may be present; this emission will, however, be localized and 
have a small surface covering factor The azimuthal averaging used will also 
tend to minimize its contribution. 

*^ T his canonical ratio of ft/zo = 5 has also been confirmed bv lKregel et alj 
12002ft who for a sample of 34 edge-on galaxies find a value of h/zo = 4.8 ± 
1.3. This is consistent with h/zo = 5 which we will continue to use. 



motivated pseudo-isothermal (ISO) halo, dominated by a cen- 
tral constant-density core. 



5.3.1. NFWhalo 
The NFW mass density distribution takes the form 

Pnfw(^) = -J (6) 

{R/R,){l+R/R,) 

where Rs is the characteristic radius of the halo and p, is re- 
lated to the density of the universe at the time of collapse of 
the dark matter halo. This mass distribution gives rise to a 
halo rotation curve 

"ln(l+cx)-cjc/(l+cjc)l'/^ 



V{R) = V200 



(7) 



x[ln(l+c)-c/(l+c)] 

where x = 7?/7?2oo- It is characterized by a concentration pa- 
rameter c = R2(K)/Rs and a velocity V2oo- These are directly re- 
lated to Rs and p, and used to parameterize the rotation curve. 
The radius R200 is the radius where the density contrast with 
respect to the critical density of the universe exceeds 200, 
roughly the virial radius (Navarro, Frenk & White 1996). The 
characteristic velocity V200 is the velocity at this radius. 

We can describe the steepness of the inner slope of the 
mass density profile with a power law p ^ r". The NFW 
models then imply a value of a = -1, the so-called "cusp". 
Note that many different pa rameterizations of CDM halos 
exist in the literature (e . g., iNavarr o. Frenk & White 19961 



1997i:lMoore et al.ll999bl:lNavarro et al..2004;.Diemand et all 



20051) . These all agree on the presence of a density cusp, but 
the steepness of the inner density profile depends slightly on 
the simulation used, with many models producing slopes that 
are steeper than the original NFW a = -l value. From an ob- 
servational point of view, however, these models all give sim- 
ilar results, and we therefore use the (in terms of inner slope, 
conservative) NFW parameterization. 

The NFW c and Vaoo par ameters are not i ndepen dent and 
are set by the cosmology. iMcGaugh et al.l ( 120071) presents 
the relation between c and V200 for two recent and represen- 
tative ACD M models. The first one is the "vanilla CDM" 
model fromlTegmark et al.'('2004'); the second one is based on 
WMAP results (Spergel et al. 2007). Both models predict a 
nearly flat distribution of c as a functio n of ¥200- Simu lations 
suggest a scatter in log(c) of ~ 0.18 ( iBullocket al.l I2OOI1) . 
The main difference between the two models is a different 
normahzation, resul ting in lower concentration values for the 
ISpergel et al.l (|2007|) model. If we consider the union of both 
models and take their respective la scatters into account, we 
expect to find values c ~ 8;!;3 for V200 = 50 km s"', to values 
c ~ for V200 = 200 km s~' . As we are considering the 
union of two models with different normalizations, this repre- 
sents a generous allowable range of c values, comparable to a 
2(7 scatter for any of the two models separately. Anticipating 
our results, we note that the c parameter indicates the amount 
of collapse a halo has undergone. Larger values of c indicate 
a larger collapse factor A value c = 1 indicates no collapse; a 
value c < 1 makes no sense within the CDM context, and we 
will refer to these as unphysical values in subsequent model- 
ing. 

5.3.2. ISO halo 

The spherical pseudo-isothermal (ISO) halo has a density 
profile 

Piso(R) = Po[l+[—) \ , (8) 



28 



de Blok et al. 



where po is the central density of the halo, and Rc the core 
radius of the halo. The corresponding dark matter rotation 
curve is given by 



V{R) = ^47rGpoRl [l " 
The asymptotic velocity of the halo, Voo, is given by 



Voo = J^nGpoRl. 



(9) 



(10) 



As the inner mass-density distribution is dominated by an 
(almost) constant-density core, it can be approximated by a 
power-law description p ^ r°' with a = Q. 

5.4. Constructing the mass models 

We use the observed rotation curve, the derived gas rotation 
curve and the rotation curve(s) of the stellar component(s) as 
input for the mass modeUng with the GIPSY task ROTMAS. 
This task subtracts the (squared) rotation curves of the bary- 
onic components (after appropriate scaling with from 
the (squared) observed rotation curve, and fits the desired dark 
matter halo model to the residuals. We use inverse-squared 
weighting of the rotation curve uncertainties in making the 
fits. For each combination of galaxy and dark matter halo 
model we derive two sets of fits. Firstly, we fix the T^^ val- 
ues of the stellar components to the values predicted by our 
empirical color-T^ ''-relations for both the diet-Salpeter and 
the Kroupa IMFs and determine the best-fitting NFW and ISO 
halo parameters. We will refer to these fits as the "fixed T^^" 
models. Secondly, we derive a set of fits where T^^ is left 
as a free parameter That is, we ignore any prior knowledge 
derived from population synthesis models, and let the fitting 
program choose its best value. These will be referred to as the 
"free T^ ''" models. 

6. RESULTS FOR INDIVIDUAL GALAXIES 

Below we describe the mass models for each galaxy along 
with a motivation for specific choices made for individual 
galaxies. The results for IC 2574 and NGC 2366 are described 
in Oh et ^1- (2008)). The T^^ values and halo model param- 
eters (including those for IC 2574 and NGC 2366) are given 
in Table |3]for the fixed T^^ models with a diet-Salpeter IMF, 
in Table |4] for a fixed Tj-^ with a Kroupa IMF, in Table |5] for 
the free T^^ ISO models, and in Table |6] for the free 
NFW models. The free models do not depend on the 
IMF. The 3.6 pm surface brightness and T^^ profiles, as well 
as the fixed and free T^^ models assuming NFW and ISO 
models are shown in Figs. 122^561 The relations of with 
(J-K) shown in these Figures are valid for a diet-Salpeter 
IMF. Values for assuming a Kroupa IMF can be obtained 
by shifting the relations down by 0.15 dex. In the interests of 
space, we only show the diet-Salpeter mass models, as they 
are generally similar to the ones derived assuming a Kroupa 
IMF. Noteworthy exceptions are discussed below in the sec- 
tions for individual galaxies. Readers more interested in the 
outcomes for the sample as a whole can at this point skip to 
Sect. 7. In the following discussion all T^ *' values assume a 
diet-Salpeter IMF, unless noted otherwise. 

6.1. NGC 925 

The surface brightness profiles of NGC 925 are shown in 
Fig. |22] The 2MASS J, H and K profiles can be traced out to 



~ 250"; the 3.6 pm profile can be traced over the entire extent 
of the Hi disk. NGC 925 shows no evidence for a bright 
central component. There is also no evidence for a strong 
color gradient, and we assume a constant T^^ = 0.65 for the 
disk. 

The rotation curves are shown in Fig. |23] The predicted 
T^-^ value make NGC 925 slightly exceed maximum disk. 
The ISO model can to some degree accommodate this, but 
the NFW model less so. This does not change when adopting 
the Kroupa T^^ values. Fits with T^^ free prefer lower values 
for T^'^. For the NFW case a negative T^^ is preferred; we 
have assumed a fixed T^'* = in this case. For both fixed and 
free Tj^ the NFW model yields parameters that are inconsis- 
tent with the CDM expectations (cf. Sect. 5.3.1), resulting in 
V200 values many times larger than the maximum rotation of 
the galaxy and c values approaching zero. This is a general 
symptom o f an NFW halo trying to fit a linearly rising rota- 
tion curve (iMcGaugh & de BloklfT998h . The ISO models in 
general fit the data better than the NFW models. 

6.2. NGC 2403 

Fig- im shows the surface brightness profiles derived from 
the 3.6 pm data, as well as the 2MASS J, H and K profiles. 
The 3.6 pva profile can be reliably traced to ^ 430" galacto- 
centric radius, the 2MASS profiles out to ~ 380". The {J-K) 
color shows a clear trend, from ~ 0.9 in the inner parts to 
^ 0.6 at the outermost measured radii. This translates into a 
trend of T^^ ~ 0.6 in the center to an outer T^-^ ~ 0.3. Be- 
yond the radius where {J-K) (and therefore T^-^) could be 
reliably determined we assume a constant value T^-* = 0.3. 
The surface brightness profile shows evidence for a second 
component in the inner parts. Based on its rather modest ap- 
pearance it is not clear though whether this constitutes a dy- 
namically separate component. We therefore consider both 
cases. For the single-component case we simply multiplied 
the surface brightness profile with the derived T^^ values. For 
the double-component case we decomposed the profile in two 
exponential profiles. The outer component is described by a 
central surface brightness =16.9 mag arcsec"^ and a scale 
length /i = 1.81 kpc. The inner component has p^-j^ = 16.7 mag 
arcsec"- and h = 0.41 kpc. We used the fitted model as a de- 
scription for the inner disk. For the outer disk we used the 
observed surface brightness distribution, with the inner disk 
model subtracted. We extended the observed outer disk pro- 
file with the corresponding fit at large radii where the disk was 
no longer directly observable. 

The {J-K) color shows a slight increase in the inner parts. 
For the inner disk we therefore assumed a constant value of 
T^ '' = 0.6, whereas for the outer, main disk we assumed the 
single-disk color distribution, the only difference being that 
we adopted a constant T^ *' = 0.5 for the innermost {R < 40") 
parts. 

Figures|25]and|26]present the fits for the single and multiple 
disk models, respectively. The NFW models (fixed and free) 
are of the same quality as the free ISO model. The fixed T^^ 
ISO models are of somewhat lower quality. NGC 2403 in 
combination with the ISO model prefers higher T^^ values 
than predicted from the infra-red color 

6.3. NGC 2841 

The surface brightness profiles of NGC 2841 are shown in 
Fig.|22l The 2MASS 7, H and K profiles can be traced out to 
- 290", the IRAC 3.6 pm profile out to - 375". All surface 
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Fig. 20. — Comparison of the NGC 7793 THINGS rotation curve with the curve derived in lCarignan & Puc"h3jl99d) . Symbols and references are indicated in 
the figure. 



brightness profiles show an extended outer exponential disk, 
as well as a more compact central component. We decompose 
the profile using a double exponential disk. For the outer disk 
we find /io = 15.7 mag arcsec"^ and h = 4.20 kpc. The in- 
ner disk has parameters /iq = 13.5 mag arcsec"^ and h = 0.72 
kpc. In subsequent modeling we use the inner disk model for 
the central component. For the outer component we use the 
total observed surface brightness profile with the inner disk 
model subtracted, extrapolated at large radii using the outer 
exponential model parameters. The surface brightness profile 
shows some evidence for central emission in excess of the ex- 
ponential model. However, the central parts of the velocity 
field cannot be modeled due to lack of H I emission. 

NGC 2841 possesses a modest color gradient, which re- 
sults in a T^^ trend from T^^ ~ 0.8 in the inner parts to 
T^-* ~ 0.5 in the outer parts. For radii beyond which the 
2MASS color could not be reliably determined, we assume 
a constant T^^ = 0.5. The color profile shows a slight red- 



dening in the innermost part, but this occurs at radii where no 
velocity information is available. 

The mass models are presented in Fig. |28] For the ISO 
model, the disk mass obtained through a free fit is very close 
to the mass derived from the population synthesis modeling 
assuming a diet-Salpeter IMF. The resulting halo has the high- 
est (well-determined) central density in the sample. The NFW 
model with fixed Tj^ is of slightly lower quality, its main 
problem is the under-prediction of the velocity at the smallest 
radii. An increase in of the inner component when T^^ 
is released compensates for this, and yields a model of equal 
quality as the ISO model. For the Kroupa fixed T^^ the ISO 
model results in an extremely compact halo, probably indicat- 
ing that T^ *" is underestimated. 

6.4. NGC 2903 

The surface brightness profiles of NGC 2903 are shown in 
Fig.|29l The 2MASS J, H and K profiles can be traced out to 
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Fig. 21 . — Comparison between the stellar disk masses derived using our method (horizontal axis) and the method described in lLerov et alj 420081) (vertical 
axis). The filled circles indicate our diet-Salpeter IMF masses, the open circles the Kroupa IMF masses. The drawn line is the line of equality, the dotted line a 
least-squares fit to the diet-Salpeter IMF masses with a slope of 0.94 and an offset of 0.44. The scatter is 0.12 dex. The dashed line shows the fit to the Kroupa 
IMF masses with a slope of 0.95, an offset of 0.39 and a scatter of 0.13 dex. 



~ 240"; the IRAC 3.6 fim profile out to - 400". The break 
in the profile around 230" is due to the sudden change in in- 
clination at this radius. The profile can be decomposed into 
two components: the inner one is modeled using an exponen- 
tial disk with /io = 13.6 mag arcsec"^ and h = 0.17 kpc. For 
the outer disk we used the observed profile with the inner disk 
profile subtracted. At radii larger than 400" we used an expo- 
nential extrapolation of the main disk. 

NGC 2903 shows a well-defined color gradient. For the 
inner disk we assumed T^'^ =1.3 which is the average value 
of TJ'^ within 20". For the outer disk we assumed Tj'^ = 0.75 
in the innermost parts and T^^ = 0.4 in the outermost parts. 
Between 20" and 150" we adopted the observed T^^ profile. 

The rotation curve models for NGC 2903 are shown in 
Fig. [30] It is immediately obvious that the inner part of 
NGC 2903 is not well fitted by the models. The models with 



T^-* fixed severely over-predict the data, independent of IMF. 
Models with T^^ as a free parameters can only be made to ap- 
proximately fit the data by setting T^-^ = for the inner disk. 
This large discrepancy is most likely due to the bar. A har- 
monic decomposition of the velocity fiel d shows significant 
non-c ircular motions in the inner 1 kpc dTrachternach et all 
l2008l) . Because of the low quality of the fits we do not list 
them in Tables [3] |5] and |6] An additional source of uncer- 
tainty could be the presence of a very dense inner molecu- 
lar disk, with signific ant amounts of associated star formation 
(iPlanesas et al.l(l997h . 

We performed an extra set of fits to those parts of the rota- 
tion curve that are unaffected by the central non-circular mo- 
tions (the flat part of the curve at R > 3.3 kpc). These are 
shown in Fig.[3T] We find that with T^ *' as a free parameter, 
the resulting disk is dynamically more significant than pre- 
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Fig. 22. — Surface brightness and T; profiles for NGC 925. Top panel: The thin full line, long dashed and short dashed curve show the 2MASS J, H and K 
surface brightness profiles, respectively. The thick full curve shows the Spitzer IRAC 3.6 /^m surface brightness distribution. None are corrected for inclination. 
Bottom panel: The thin full curve shows the observed 2MASS (J—K) distribution, with the thin dotted line indicating the average value. The thick full line shows 
the distribution of derived from (J—K) as described in the text. The thick dashed curve shows our assumed model for the Tj-* distribution of the disk. 



dieted by population synthesis models and (J-K) colors. The 
ISO model prefers a dynamically insignificant central com- 
ponent, whereas the NFW model prefers a central T^^ value 
close to the predicted one. Here also, the main disk is found 
to be more massive than predicted by the colors. The NFW 
model fits better than the ISO model but due to the presence 
of the non-circular motions, it is not clear whether this can be 
generalized to the entire rotation curve. For fixed T^^, the 
ISO model demands an extremely small core radius Rc and a 
very large central density po, in order to produce the flat halo 
rotation curve shown in Fig. [311 



6.5. NGC 2976 



The surface brightness profiles of NGC 2976 are shown in 
Fig. [32l The 2MASS J, H and K profiles can be traced out 
to ^ 210". The H I disk and rotation curve extend only out 
to 150", and no extrapolations using exponential disks were 
necessary. The galaxy shows a well-defined color gradient, 
leading to a T^^ variation from T^^ ~ 0.75 in the inner parts 
to T^-* ~ 0.3 in the outer parts. We adopted the T^^ distribu- 
tion as shown in Fig.[32]without any further modifications. 

The rotation curve fits are shown in Fig. [33] The models 
with the fixed diet-Salpeter Tj^ value slightly over-predict 
the observed rotation curve. The free-T^^, ISO model also 
prefers a model very close to maximum disk, with a value for 
T^-* close to the fixed T^ '' Kroupa value. The large value 
of Rc, seems to indicate NGC 2976 is located well within 




Fig. 23. — ISO and NFW rotation curve fits for NGC 925. Ttie top row siiows fits using tlie ISO halo. The bottom row those assuming the NFW model. The 
left-hand column shows the fits assuming fixed (predicted) values of T^^. The right-hand column shows the fits produced with a free In all panels the dots 
represent the observed curve and the thin full lines the uncertainties. The grey dotted curve shows the rotation curve of the gas; the thin dashed grey curve that of 
the stellar disk. The thick, grey full curve shows the resulting rotation curve of the halo. The thick long-short dashed curve the resulting best-fit model rotation 
curve. 



the core radius of its halo. The free-Tj-^, NFW fit yields un- 
physical halo parameters trying to approximate a solid-body 
curve. It can only do this by choosing a large value of V200 
and a small value of c. As mentioned before, c represents a 
collapse factor, and it is clear that a value c < 1 makes no 
physical sense. Similarly, it is unlikely that a fairly modest 
galaxy such as NGC 2976 is embedded in a halo with a typi- 
cal velocity V200 > 500 km s"' . 

6.6. NGC 3031 

The surface brightness profiles of NGC 303 1 are shown in 
Fig. |34] The 2MASS /, H and K profiles can be traced out to 
~ 750"; the 3.6 fim profile out to ~ 850". The profiles show 
clear evidence for a central component. The outer disk shows 



a gradual flattening of the surface brightness profile towards 
larger radii. We modeled the central component using an ex- 
ponential disk with parameters /io = 12.2 mag arcsec"^ and 
h = 0.25 kpc. To describe the outer disk we simply subtracted 
the central exponential disk from the total light profile. The 
slight 0.2 mag) residuals seen around R ^ 100" do not im- 
pact on any subsequent modeling. The (J-K) color profile 
shows no significant color gradient, except for a central red- 
der component. We assume a constant T^-^ = 0.8 for the disk 
and T^-* =1.0 for the central component. 

The rotation curve mass models are shown in Fig. [35] Al- 
ready at first glance, it is clear that the model fits deviate sig- 
nificantly from the observed rotation curve. Although the gen- 
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Fig. 24. — Surface brightness and T^^ profiles for NGC 2403. Lines and symbols ai'e as in Fig.|22] In addition, the thin gray lines in the top panel show the 
decomposition of the profile in multiple exponential disk models. The long-dashed grey line in the bottom panel shows the assumed value of T^'^ for the inner 
disk. 



eral trends in the curve are described well, smaller scale devi- 
ations clearly indicate the presence of significant non-circular 
motions. Given the presence of the prominent, grand-design 
spiral arms (causing the bump at R ^ 7.5 kpc), and MSl's 
location in an interacting system, this should not come as a 
surprise. The fit with fixed T^^ yields values that are only 
slightly inferior to the free T^^ fits. 

6.7. NGC 3198 

The surface brightness profiles of NGC 3198 are shown in 
Fig. |36] The 2MASS /, H and K profiles can be traced out to 
~ 240"; the 3.6 /im profile out to ^ 250". The profiles show 
a central component, as well as a "shoulder" in the surface 
brightness around R ~ 50". We initially model the system us- 
ing two components. For the innermost component we use an 



exponential disk with parameters yuo =15.4 mag arcsec"^ and 
h = 0.56 kpc. The profile of the main disk is derived by sub- 
tracting the inner component model from the total luminosity 
profile and is extended beyond R ~ 240" by an exponential 
disk model with /iq = 16.2 mag arcsec"- and h = 3.06 kpc. 

The most striking feature in the (J-K) color distribution is 
the sudden reddening in the innermost couple of arc-seconds. 
At face-value this would imply a steep increase in T^'^ in 
the innermost parts. This might indicate an additional (third) 
compact, central component. Alternatively, it could indi- 
cate the presence of a centrally concentrated pronounced dust 
component. Whatever the cause, the steep increase in color 
only affects the very innermost point of the rotation curve and 
we will not model it. Adding a third component would in- 
crease the number of possible parameter combinations, but 
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Fig. 25. — ISO and NFW rotation curve fits for the single-disk model of NGC 2403. Symbols and lines as in Fig. 1231 



without obvious benefit. We therefore adopt the T^^ profile 
as observed, but assume a constant T^ *" = 0.8 for R > 125", 
as well as a constant T^^ = 0.7 for R < 10", as indicated in 
Fig. [361 

The rotation curve fits are shown in Fig. [37] In the case 
of fixed TJ^, the inner component clearly overestimates the 
rotation velocity, indicating that it is likely to be dynamically 
less important than suggested by the models (assuming the 
rotation curve is not affected; see below). With T^^ as a free 
parameter the outer disk ends up being less massive than pre- 
dicted for both models. The central component is also fitted 
with a much smaller Tj^. For the ISO model we find an in- 
ner component ~ 10 times lighter than predicted, whereas the 
NFW model results in a negative, hence unphysical, value for 
T|-^. We therefore fixed its value to = 0, with little im- 
pact on the results for the outer disk. We checked whether 



modeling the inner component as a spherical bulge instead of 
an exponential disk changed the results, but found this had a 
negligible effect on the outcome. Because of the apparent in- 
significance of the inner component, we also investigated if a 
single-disk model could explain the data. The results are also 
given in Fig. [38] It is clear that independent of the choice for 
halo model, the inner radii cannot be well fitted. 

The small value of for the inner disk makes the free 
T^-^ fits less satisfactory. The low T^^ values in the inner 
parts seems to contradict the sudden reddening in the inner 
10", suggesting a higher Tl'^ value. A close inspection of the 
IRAC 3.6 fim image indicates the presence of a short bar at 
a small angle with the minor axis, which could mean that the 
rotation velocities are affected by non-circular motions, thus 
giving incorrect values for T^^. However, a harmonic decom- 
position of the velocity field shows only small non-circular 
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Fig. 26. — ISO and NFW rotation curve fits for tlie multiple-disic model of NGC 2403. Symbols and lines as in Fig. 1231 In addition, the long-dashed grey 
curve represents the rotation curve of the inner disk. 



mo tions (of order a few kms~ ') in the inner parts of the disk; 

see lTrachternach et al.l (l2008h . 

In an independent analysis, iBottema et al.l (|2002|) also note 
the presence of a possible bar in a K' image of NGC 3198. 
They argue that having a bar roughly parallel to the minor 
axis would increase the apparent rota tion velocities in the in- 
ner region ( iTeuben & Sanderslll985l) . If this is indeed what 
is happening in NGC 3198, then the true rotation velocity 
in the inner parts must be even lower than we now observe, 
thus making the inner disk even less massive. Clearly there is 
scope here for a more extensive study of the dynamics in the 
inner part of NGC 3198. 

6.8. NGC 3521 

The surface brightness profiles of NGC 3521 are shown in 
Fig. |39] The 2MASS 7, H and K profiles can be traced out 



to ^ 250"; the 3.6 profile out to the edge of the H I disk. 
The photometric analysis is made difficult by the faint stellar 
halo surrounding the main disk. It is significantly extended 
along the minor axis, which, through the ellipse integration 
used, results in a flattening of the surface brightness profile at 
larger radii. For this reason results from the stellar decompo- 
sition are more uncertain than for other galaxies in the current 
sample. 

We treat the stellar disk as a single component. A clear 
color gradient is present. The value for changes from 
~ 1 in the center to ~ 0.6 in the outer parts. Beyond the last 
reliable 2MASS radius we assume a constant T^^ = 0.6. The 
rotation curve models are presented in Fia.|40] Neither model 
fits perfectly, but both give reasonable T^ *" ratios. Assuming 
the fixed T^-^ values, it follows that the dynamics of NGC 
3521 is dominated by its stars. 
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Fig. 27. — Surface brightness and T^^ profiles for NGC 2841. Lines and symbols as in Fig. 1241 



6.9. NGC 3621 

The surface brightness profiles of NGC 3621 are shown in 
Fig. EI] The 2MASS 7, H and K profiles can be traced out 
to 320"; the 3.6 profile out to 400". The latter pro- 
file shows slight evidence for a change in scale length around 
R ^ 200". However, neither the lack of a clear break in the 
color gradient, nor the morphology, nor the (small) deviation 
from a single exponential disk are enough justification for a 
two disk model. We extended the disk to larger radii using 
an exponential fit to radii 200" < R < 400" and with param- 
eters /io = 16.6 mag arcsec"^ and h = 2.61 kpc. A small color 
gradient is visible in the 2MASS data which results in a gra- 
dient in Tl^ from T^^ ~ 0.8 in the inner parts to T^^ ~ 0.4 
in the outer parts. Due to the increased uncertainty in the 
2MASS color measurement in the outermost parts, we as- 
sumed Tl^ = 0.45 for R > 150", as indicated in Fig. ED 



The rotation curve fits are presented in Fig. |42l It is inter- 
esting how well the best (free) fit and the fixed diet-Salpeter 
fit for the ISO halo agree with each other Similarly, the free 
NFW fit yields a value of T^-^ very close to the Kroupa fixed 
T^ *" value, a clear illustration of how the NFW model prefers 
lighter disks in order to accommodate its steeper halo mass 
density profile. The free ISO and NFW fits are of identical 
quality, and on the basis of these data it is difficult to single 
out a favorite model for NGC 3621. 

6.10. NGC 4736 

The surface brightness profiles of NGC 4736 are shown in 
Fig. EH The 2MASS J, H and K profiles can be ti-aced out to 
~ 250", whereas the 3.6 /im profile can be traced out to the 
edge of the H I disk at ~ 400". 

The stellar component of NGC 4736 has a complex struc- 
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Fig. 28.— ISO and NFW rotation curve fits for NGC 2841. Lines and symbols as in Fig.l26l 



ture. A prominent stellar ring can be seen in various wave- 
length at a radius R ^ 40". The ring, visible in the 3.6 /im 
image, is also prominently visible in the 20-cm radio contin- 
uum. The presence of this fairly strong continuum suggests 
active star formation, and may mean that at this particular ra- 
dius the assumption that the 3.6 /im emission is proportional 
to the stellar mass is not valid. For this reason we will not 
take this feature into account in our subsequent analysis of the 
surface brightness profile. The IRAC 3.6 /im profile flattens 
out at 250" < < 400", caused by the faint stellar ring sur- 
rounding the main body of NGC 4736. At larger radii (beyond 
the edge of the H I disk) the profile starts decreasing again. 
Also visible is a significant steepening in the inner parts. We 
will treat this as evidence for a separate central component 
with parameters /io = 11-8 mag arcsec"^ and h = 0.26 kpc. To 
model the outer disk we fitted an exponential disk to the ra- 



dial range 80" <R < 240" and find parameters /io = 15.9 mag 
arcsec"- and li = 1 .99 kpc. We used the extrapolated fit for 
radii R < 80" and the observed radial surface brightness for 
R > 80". 

NGC 4736 exhibits a small color gradient in the outer parts 
with a steeper inner reddening. We have adopted a constant 
value of TJ-^ = 0.6 for the outer disk beyond the last reliable 
2MASS (J-K) value, and have additionally assumed a con- 
stant TJ^ = 0.7 for the innermost 20". For the central disk we 
assume a constant T^-^ = 0.9 which is the average value for 
R > 20". 

The rotation curves are shown in Fig.|44l It is immediately 
obvious that NGC 4736 must be close to maximum disk. The 
predicted T^^ value for the inner disk yields a curve that over- 
estimates the rotation velocity (not shown here). We there- 
fore left T^^ of the inner disk as a free parameter after fixing 
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Fig. 29. — Surface brightness and T^^ profiles for NGC 2903. Lines and symbols as in Fig. 1241 



the T^^ value of the main disk. This resulted in an inner 
disk mass about a factor of ^ 2 less than predicted (and is 
the model shown in Fig. l44l i. Whilst the discrepancies are 
slighly less for the Kroupa T^-^ values, the resulting models 
still overpredict the rotation velocity. Because of these uncer- 
tainties and because of the large non-circular motions in this 
galaxy (jTrac hternach et al. 20081), it is virtually impossible to 
say anything definite about the distribution of dark matter in 
this galaxy based purely on rotation curve arguments. 

6.11. DD0154 

DDO 154 is barely detected in the 2MASS images and we 
cannot therefore use the same (J-K) color method as with 
the other galaxies. We therefore use the (B-V) a nd (B-R) 
colors as proxies, as described in lOh et al.l (l2008h. For th e 
color we use the values given in lCarignan & BeauUeul ( ll989l) : 



(B-V) = 0.37 ± 0.06 and (B-R) = 0.64 ± 0.07. No significant 
color gradient is detected, and we therefore assume a constant 
colo r as function of r adius in our models. Using the equations 
from lOh et all (l2008h we find values T^ '' = 0.34 ± 0.04 using 
the (B-V) colors and Tj^ = 0.31 ±0.02 using the (B-R) 
colors. We adopt the weighted average Tj^ = 0.32 as our 
best estimate. The IRAC 3.6 /im profile can be described by 
a single exponential disk with yUo = 20.8 mag arcsec"^ and 
h = 0.72 kpc. 

The rotation curves are shown in Fig. |46l The ISO model 
consistently fits better than the NFW model. The latter over- 
predicts the velocity in the inner parts. Note the very large 
values for T^-^ when it is left as a free parameter These will 
be discussed in more detail in Sect. 7. Because of the dark 
matter dominance in DDO 154 the derived halo parameters 
are virtually independent from the assumed IMF. 
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Fig. 30. — ISO and NFW rotation curve fits for NGC 2903 for the entire radial range. Lines and symbols as in Fig. 1261 



6.12. NGC 5055 

The surface brightness profiles of NGC 5055 are shown in 
Fig. |47] The 2MASS 7, H and K profiles can be traced out 
to ~ 350"; the 3.6 jiva profile out to 450". There is clear ev- 
idence for a compact central component which we model as 
an exponential disk with parameters /io = 13.4 mag arcsec"^ 
and h = 0.35 kpc. For the outer disk, we simply used the ob- 
served profile with the inner component subtracted, as indi- 
cated in Fig. |47l Small residuals are present, notably around 
R ^ 50", but these are of minor importance. An exponen- 
tial disk fit to radii 350" <R< 450" was used to extend the 
outer disk profile. NGC 5055 has a significantly redder cen- 
ter, with a well-defined color gradient at larger radii. For the 
central component we assume a constant T^'^ = 1.3, whereas 
for the outer disk we find values from an inner T^-^ ~ 1 to an 



outer T^ '' ~ 0.5. Beyond the radii where the 2MASS colors 
could be reliably determined we assume a constant value of 
T^^ = 0.5. In the very inner parts we assumed an extrapo- 
lated = 1 .0 for the outer disk. 

The rotation curve fits are presented in Fig.|48l In the case 
of fixed T^-^ the inner disk over-predicts the curve by a large 
factor (not shown here), and we kept it as free parameter, 
yielding values of T^-^ ~ 0.2 for the inner disk. The NFW 
model seems to have particular difficulties fitting this galaxy. 
With TJ^ as a free parameter, the fit demands a negative value 
for the outer disk. In the fits presented here we have therefore 
set T^^ = for this component. 

6.13. NGC 6946 

The surface brightness profiles of NGC 6946 are shown in 
Fig.EH The 2MASS 7, H and K profiles as well as the 3.6 ^m 
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profile can all be traced out to ^ 360". The profiles show clear 
evidence for a compact central component. We modeled this 
component as an exponential disk with parameters /io = 12.8 
mag arcsec"^ and li = 0.15 kpc. For the outer disk we used 
the observed profile with the inner component subtracted. For 
larger radii we extended the profile with an exponential fit 
with parameters /io = 16.3 mag arcsec"^ and h = 2.97 kpc. 
There is evidence for a slight excess of light in the "shoulder" 
of the profile at R ^ 50". This small excess does, however, 
have negligible effect on the results of the rotation curve fit, 
and we therefore use a two-component model for the stellar 
disk. 

The outer disk of NGC 6946 shows no strong evidence for 
a (J-K) color gradient and we use a constant Tj-^ = 0.64. 
Within R ^ 30", the change in color is, however, very pro- 
nounced, resulting in very red central components [{J-K) > 
1.4 at the innermost point]. Taking these colors at face-value 
would imply T^^ > 2.0. Such col ors are only fou nd for ex- 
treme star formation histories (cf. lOh et al.ll2008h . and it is 
therefore likely that other effects contribute to this extreme 
color For the purposes of the rotation curve analysis we there- 
fore assume T^-^ =1.0 for the inner disk (which is also the 
value one gets when averaging over the entire radial extent of 
the stellar component). 

The halo model fits are presented in Fig. |50l They show 
that the inner disk component prefers a value for that is 
lower than the (J-K) color would suggest. The values derived 
for the outer disk are reasonably close to the predicted values. 
The NFW and ISO models all produce very similar quality 
fits. 

6.14. NGC 7331 

The surface brightness profiles of NGC 7331 are shown in 
Fig. im The 2MASS J, H and K profiles can be traced out 
to ~ 265"; the 3.6 /xm profile can be traced out to 350", or 
the entire extent of the H I disk. The profiles show clear ev- 
idence for a compact central component. We modeled this 
component as an exponential disk with parameters /io = 12.0 
mag arcsec"^ and h = 0.32 kpc. For the outer disk we sim- 
ply used the observed profile with the inner component sub- 
tracted, as indicated in Fig. [51] NGC 7331 has a well-defined 
color gradient, and shows some of the largest T^^ values, as 
well as one of the steepest T\'^ gradients in the entire sample. 
Beyond the radii where the 2MASS colors could be reliably 
determined we assume a constant value of T^^ = 0.7. 

These high T^'^ values lead to incompatible mass models: 
using the predicted T^^ values results in disks that are too 
massive. This is shown in Fig. |52] clearly suggesting that 
for this galaxy the color-based values are not correct (for that 
reason we do not list the model parameters in Tables [3] |5] and 
|6]). Leaving T^-^ free in the fit, the data clearly prefer lower 
values. A possible explanation could be the p resence of a 
strong dust ring in the inner disk of NGC 7331 dRegan et alj 
120041) . The radius at which this ring is found corresponds 
with the radius where the highest T^^ values are found. This 
would suggest that the very inner Tj^ maximum and steep 
drop are associated with the central component, whereas the 
subsequent steep rise and gradual drop are associated with 
the dust ring, and therefore do not reflect stellar population 
changes. 

We therefore also investigate a model where the two com- 
ponents each have a radially constant Tj^. Using Fig.|5T] we 
find that T^^ = 1.0 and Tj '' = 0.7 are good estimates for the 



inner and outer disks, respectively. These fits are presented in 
Fig. |53] We then find that both ISO and NFW models result 
in better fits. NGC 7331 is therefore the one galaxy in our 
sample where the {J—K)-Tl:^ relation clearly fails. Note that 
for the fixed T^^ case, the Kroupa models yield much better 
fits than the diet-Salpeter ones. 

6.15. NGC 7793 

The surface brightness profiles of NGC 7793 are shown in 
Fig. m The 2MASS /, H and K profiles can be traced out 
to ~ 325"; the 3.6 /im profile can be traced out to 400" - 
the entire extent of the Hi disk. The 2MASS profiles show 
a clear and sudden increase in surface brightness in the very 
inner part. This is the signatu re of the nuclear star cluster in 
NGC 7793 jBoker et al.li2002l) . Its photometric and dynami- 
cal importance has already become negligible at the radius of 
the innermost point of the rotation curve. The extent of the 
IRAC 3.6 /im profile equals that of the Hi data, and no ex- 
ponential extrapolations were necessary. The {J-K) profile 
shows a small but consistent color gradient, translating in a 
gradient in T^^ from ~ 0.5 in the inner parts, to ~ 0.25 at 
the outermost radius where the 2MASS colors were deemed 
reliable. Beyond this radius we assume a constant T^-* = 0.25 
as indicated in Fig.|54l 

The rotation curve mass models are presented in Fig. |55] 
For the fixed T^-^ model both ISO and NFW models give 
roughly equal results: neither is fully able to describe the ob- 
served rotation curve (independent of IMF). The "bump" in 
the rotation curve within R = 2 kpc, and its resemblance with 
the shape of the stellar disk rotation curve suggests that the in- 
ner part of NGC 7793 could be close to maximum disk. The 
ISO model with T^-* as a free parameter does indeed prefer 
T^^ values that are slightly higher than the fixed T^-^ diet- 
Salpeter value. For the NFW model, releasing T^^ has forced 
its value down to unrealistic values in order to accommodate 
the more concentrated dark halo. 

As remarked earlier (Sect. 14.171 ). the possible observed de- 
cline in the outer rotation curve is uncertain. We tested 
whether this decline could cause the low quality of the fits by 
repeating the analysis with radii R > 4.7 kpc excluded. This 
analysis is shown in Fig. |56l For the fixed case the sit- 
uation does not change much. The biggest improvement is 
achieved for the ISO model with a free T^ ''. The preferred 
value is close to maximum disk, with a significantly better 
overall fit. 

7. DISCUSSION 

We now put the results derived for individual galaxies in a 
broader context by looking at the sample as a whole. Sect. 7. 1 
deals with some of the general properties of the rotation 
curves, and looks at correlations with luminosity, as well as 
the mass fractions within the galaxies. Sect. 7.2 deals with 
the values for Tj'' and compares them with other measures 
for the stellar mass-to-light ratio. Sections 7.3 and 7.4 deal 
with the halo properties. In Sect. 7.3 we discuss the halo rota- 
tion curves, and in Sect. 7.4 the parameters of the halo models 
are compared with predictions based on cosmological models. 

7.1. General properties 

THINGS contains galaxies spanning a large range of lu- 
minosities and Hubble types. Consequently, there is a large 
range in shape, amplitude and extent of the rotation curves 
within our sample. To illustrate this diversity, we plot in 
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TABLE 3 

Mass Models with fixed T^-^and diet-Salpeter IMF 



ISO halo NFWhalo 



Name logMf (T3-6> logMf (T3-|) Rc po x? c V200 X? 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


(9) 


(10) 




(11) 


NGC 925 


10.01 


0.65 






16.86 ± 7.47 


3.4 ± 


0.5 


2.15 


<0.1 


>500 




2.81 


NGC 2366 


8.41 


0.33 






1.36 ± 0.07 


34.8 ± 


2.4 


0.17 


<0.1 


>500 




0.98 


NGC 2403 (1 comp) 


9.71 


0.41 






2.09 ± 0.05 


81.2 ± 


3.6 


0.88 


9.9 ± 0.2 


109.5 ± 


1.0 


0.55 


NGC 2403 (2 comp) 


9.67 


0.39 


8.63 


0.60 


2.14 ± 0.05 


77.9 ± 


3.3 


0.79 


9.8 ± 0.2 


110.2 ± 


1.0 


0.56 


NGC 2841 


1 1.04 


0.74 


10.40 


0.84 


2.03 ± 0.05 


298.7 ± 


14.9 


0.27 


16.1 ± 0.2 


183.2 ± 


1.2 


0.42 


NGC 2903 (outer) 


10.15 


0.61 


9.33 


1.30 


<0.01 ± 0.09 


>1000 




0.63 


30.9 ± 0.6 


112.9 ± 


0.6 


0.36 


NGC 2976 


9.25 


0.55 






>1000 


11.7 ± 


2.1 


1.76 


<0.1 


>500 




2.78 


NGC 3031 


10.84 


0.80 


10.11 


1.00 


5.25 ± 1.36 


14.8 ± 


4.2 


3.93 


3.0 ± 2.9 


190.9 ± 


161.1 


4.36 


NGC 3198 (1 comp) 


10.40 


0.80 






3.22 ± 0.16 


33.5 ± 


3.0 


0.84 


7.5 ± 0.4 


112.4 ± 


2.1 


1.37 


NGC 3198 (2 comp) 


10.45 


0.80 


9.46 


0.73 


4.97 ± 0.41 


14.4 ± 


2.0 


2.15 


5.1 ± 0.5 


122.7 ± 


4.9 


2.88 


IC 2574 


9.02 


0.44 






7.23 ± 0.30 


4.1 ± 


0.1 


0.17 


<0.1 


>500 




1.81 


NGC 3521 


11.09 


0.73 






39.4 ± 90.4 


1.3 ± 


1.2 


8.07 


<0.1 


403.2 ± 


123.2 


8.52 


NGC 3621 


10.29 


0.59 






5.54 ± 0.16 


14.4 ± 


0.6 


0.62 


3.7 ± 0.2 


165.5 ± 


5.9 


0.81 


NGC 4736 


10.27 


0.63 


9.59" 


0.33" 


1.44 ± 1.57 


22.4 ± 


41.5 


1.52 


11.4 ±9.8 


35.2 ± 


0.3 


1.51 


DDO 154 


7.42 


0.32 






1.34 ±0.06 


27.6 ± 


1.6 


0.44 


4.4 ± 0.4 


58.7 ± 


4.3 


0.82 


NGC 5055 


11.09 


0.79 


9.32'' 


0.11'' 


45.63 ± 0.24 


0.9 ± 


0.2 


8.13 


<0.1 


450.1 ± 


32.4 


10.31 


NGC 6946 


10.77 


0.64 


9.58 


1.00 


20.58 ± 3.77 


5.4 ± 


0.4 


1.45 


<0.1 


>500 




2.59 


NGC 7331 (const)' 


11.22 


0.70 


10.24 


1.00 


>1000 


1.6 ± 


0.1 


2.93 


<0.1 


>500 




4.08 


NGC 7793 


9.44 


0.31 






1.93 ±0.15 


77.2 ± 


7.7 


2.97 


5.8 ± 1.4 


156.6 ± 


39.1 


4.17 


NGC 7793 (rising) 


9.44 


0.31 






3.52 ± 0.51 


50.9 ± 


4.8 


2.26 


<0.1 


>500 




3.67 



Note. — (2): Logarithm of the predicted stellar mass of the disk (Mq ). (3): Average stellar mass-to-light ratio in the 3.6 jum band (Mq /Lq ). (4): Logarithm of 
the predicted stellar mass of the bulge (Mq ). (5): Average stellar mass-to-light ratio of the bulge component in the 3.6 band (Mq /Lq ). (6) Core radius Rc and 
associated uncertainty (kpc). (7) Central density po and associated uncertainty (IO^'Mq pc"'). (8) Reduced xl- (9) NFW c parameter and associated uncertainty. (10) 
NFW parameter V200 and associated uncertainty (km s"'). (11) Reduced xl- 
" JWbuige is a free parameter. Predicted model values are Tj^ = 0.9 and logMj = 10.02. 

A^buige is a free parameter. Predicted model values are T''^ = 1.3 and logAfJ = 10.40. 

The model with the color gradient included severely over-predicts the data and is not Ustedhere. Its parameters are = 0.96 and logM° = 11.36 for the main disk 
and T^-^ = 1.08 and logM* = 10.28 for the inner disk. 



TABLE 4 

Mass Models with fixed TJ ^and Kroupa IMF 



ISO halo NFWhalo 



Name 
(1) 


logMf 
(2) 


(3) 


logMf 
(4) 


(5) 


Rc 
(6) 


PO 

(7) 




(8) 


c 

(9) 


^200 

(10) 


(11) 


NGC 925 


9.86 


0.47 






9.67 ± 1.27 


5.9 ± 


0.5 


1.14 


<0.1 


>500 


2.17 


NGC 2366 


8.26 


0.23 






1.32 ± 0.07 


37.3 ± 


2.4 


0.16 


<0.1 


>500 


1.01 


NGC 2403 (1 comp) 


9.56 


0.29 






1.49 ± 0.05 


152.8 ± 


7.5 


1.04 


12.4 ± 0.2 


101.7 ± 0.7 


0.57 


NGC 2403 (2 comp) 


9.52 


0.26 


8.48 


0.43 


1.52 ± 0.04 


145.8 ± 


6.9 


0.97 


12.3 ± 0.2 


102.2 ± 0.7 


0.57 


NGC 2841 


10.88 


0.53 


10.25 


0.60 


0.63 ± 0.04 


3215.3 ± 


371.8 


0.22 


18.9 ± 0.4 


181.4 ± 1.0 


0.23 


NGC 2903 (outer) 


10.00 


0.43 


9.18 


0.92 


<0.01 ± 0.09 


>1000 




1.14 


35.5 ± 0.7 


111.8 ± 0.6 


0.41 


NGC 2976 


9.10 


0.39 






5.09 ± 2.54 


35.5 ± 


3.1 


0.50 


<0.1 


>500 


1.90 


NGC 3031 


10.69 


0.57 


9.96 


0.71 


0.78 ± 0.19 


754.2 ± 


323.6 


3.88 


26.4 ± 2.5 


94.6 ± 3.9 


3.61 


NGC 3198(1 comp) 


10.25 


0.57 






2.72 ± 0.13 


46.9 ± 


4.0 


0.80 


8.7 ± 0.4 


109.7 ± 1.7 


1.30 


NGC 3198 (2 comp) 


10.30 


0.57 


9.31 


0.52 


2.82 ± 0.19 


44.0 ± 


5.1 


1.41 


8.5 ± 0.5 


110.4 ± 2.2 


2.06 


IC 2574 


8.87 


0.31 






6.18 ± 0.21 


5.0 ± 


0.1 


0.17 


<0.1 


>500 


1.73 


NGC 3521 


10.94 


0.52 






2.50 ± 0.66 


73.0 ± 


30.6 


4.75 


8.9 ± 2.0 


128.4 ± 16.4 


5.55 


NGC 3621 


10.14 


0.42 






2.77 ± 0.10 


48.9 ± 2.8 


1.09 


7.8 ± 0.2 


120.2 ± 1.4 


0.55 


NGC 4736 


10.12 


0.44 


9.56" 


0.31" 


<0.01 ± 0.65 


>1000 




1.73 


63.5 ± 24.2 


42.4 ± 1.7 


1.41 


DDO 154 


7.27 


0.23 






1 .32 ± 0.06 


28.5 ± 


1.7 


0.44 


4.5 ± 0.4 


58.0 ± 4.1 


0.83 


NGC 5055 


10.94 


0.56 


9.81'' 


0.34'' 


11.73 ± 0.71 


4.8 ± 


0.4 


1.03 


2.1 ± 0.4 


217.8 ±21.2 


1.45 


NGC 6946 


10.62 


0.45 


9.43 


0.71 


3.62 ± 0.16 


45.7 ± 


3.0 


0.98 


6.2 ± 0.5 


183.8 ± 11.1 


1.03 


NGC 7331 (const)'^ 


11.07 


0.50 


10.09 


0.71 


5.40 ± 0.31 


24.4 ± 


2.1 


0.31 


4.9 ± 0.4 


200.0 ± 10.7 


0.24 


NGC 7793 


9.29 


0.22 






1.46 ± 0.10 


126.0 ± 


12.2 


3.56 


9.1 ± LI 


114.1 ± 12.2 


3.74 


NGC 7793 (rising) 


9.29 


0.22 






1.98 ± 0.22 


95.7 ± 


10.7 


4.08 


<0.1 


>500 


2.06 


Note. — (2): Logarithm of the predicted stellar mass of the disk (Mr:)). (3): Average stellar mass-to-light ratio in the 3.6 band (Mq/Lq). (4): Logarithm of 
the predicted stellar mass of the bulge (Mq ). (5): Average stellar mass-to-light ratio of the bulge component in the 3.6 pm band (Mq /Lq). (6) Core radius Rc and 



associated uncertainty (kpc). (7) Central density po and associated uncertainty (IQT^Mq pc'^). (8) Reduced x^. (9) NFW c parameter and associated uncertainty. (10) 
NFW parameter ^200 ^ind associated unceitainty (km s~'). (1 1) Reduced x^- 
^ Mbuiui; is a free parameter Predicted model values are T^'^ = 0.6 and logMf = 9.87. 
^ Mbuige is a free parameter. Predicted model values are T^^ = 0.8 and logM^ = 10.25. 

^ The model with the color gradient included severely over-predicts the data and is not listed here. Its parameters are T^"^ = 0.61 and logM^ = 11.21 for the main disk and 
T^^ = 0.69 and logA/f = 10.13 for the inner disk. 
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TABLE 5 
ISO Mass Models with free TJ * 
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(1) 


M" 
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(3) 


f 
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NGC 303 1 


10.93 


0.99 


1.24 ± 0.54 


9.54 


0.22 


0.27 ± 0.56 


4.13 ± 23.4 


12.2 ± 


124.6 


3.26 


NGC3198(lcomp) 


10.34 


0.70 


0.88 ± 0.07 








2.71 ± 0.33 


47.5 ± 


11.4 


0.81 


NGC 3198 (2 comp) 


10.30 


0.57 


0.71 ± 0.03 


8.36 


0.06 


0.08 ± 0.04 


1.86 ±0.12 


97.1 ± 


12.2 


0.36 


IC 2574 


8.37 


0.10 


0.23 ±0.18 








4.99 ± 0.34 


6.7 ± 


0.7 


0.16 


NGC 3521 


10.77 


0.34 


0.47 ±0.13 








1.32 ± 0.76 


370.2 ± 


451.1 


4.04 


NGC 3621 


10.30 


0.60 


1.02 ± 0.02 








5.88 ± 0.32 


13.0 ± 


1.1 


0.61 


NGC 4736 


10.23 


0.58 


0.92 ± 0.06 


9.56 


0.32 


0.35 ± 0.09 


0.04 ± 1.23 


>500 




1.53 


DDO 154 


8.45 


3.46 


10.82 ± 1.05 








2.69 ± 0.24 


9.0 ± 


1.1 


0.28 


NGC 5055 


10.88 


0.24 


0.62 ± 0.02 


9.69 


0.40 


0.31 ±0.01 


7.15 ± 1.12 


11.1 ± 


3.1 


0.87 


NGC 6946 


10.58 


0.41 


0.64 ±0.15 


9.50 


0.83 


0.83 ± 0.06 


3.32 ± 1.23 


55.6 ± 


42.3 


0.96 


NGC 7331 (grad) 


10.88 


0.39 


0.33 ± 0.04 


10.22 


0.95 


0.88 ± 0.08 


3.88 ± 0.98 


58.0 ± 


27.8 


0.26 


NGC 7331 (const) 


11.10 


0.32 


0.76 ± 0.03 


10.14 


0.79 


0.79 ± 0.09 


9.35 ± 1.43 


10.1 ± 


2.4 


0.26 


NGC 7793 


9.53 


0.38 


1.22 ±0.14 








2.40 ± 0.48 


53.3 ± 


15.7 


2.83 


NGC 7793 (rising) 


9.57 


0.42 


1.35 ± 0.04 








oo 


25.7 ± 


1.0 


1.01 



Note. — (2): Logarithm of the predicted stellar mass of the disk (Mq). (3): Average stellar mass-to-light ratio in the 3.6 fim band (Mq /Lq). 
(4): Scaling factor / with respect to the predicted diet-Salpeter T_^^ value listed in Table 3. The factor with respect to the Kroupa value listed 
in Table 4 can be obtained by multiplying with a factor of 1.41. (5): Logarithm of the predicted stellar mass of the bulge (Mq). (6): Average stellar 
mass-to-light ratio in the 3.6 /im band (Mq/Lq). (7): Scaling factor/* with respect to the predicted diet-Salpeter T''^ value hsted in Table 3. The 
factor with respect to the Kroupa T^ ' value listed in Table 4 can be obtained by multiplying with a factor of 1.41. (8) Core radius Rc and associated 
uncertainty (kpc). (9) Central density po and associated uncertainty (W'^Mq pc~'). (10) Reduced xl- 
" The stellar mass-to-light ratio T^'^ was fixed to T^'^ = 0. 



TABLE 6 

NFW Mass Models with free TJ * 



Name 








Mf 






c 




^200 






(1) 


(2) 


O) 


(4) 


O) 


(6) 


(7) 


(S) 




('» 




(10) 


NGC 925 




0.00» 


0.00» 








<0.1 




>500 




1.11 


NGC 2366 




0.00^ 


0.00=' 








<0.1 




>500 




1.44 


NGC 2403 (1 comp) 


9.69 


0.39 


0.96 ± 0.05 








10.2 ± 


0.5 


108.3 ± 


1.9 


0.55 


NGC 2403 (2 comp) 


9.59 


0.33 


0.83 ± 0.08 


8.62 


0.60 


0.99 ± 0.09 


10.9 ± 


0.6 


106.1 ± 


1.9 


0.55 


NGC 2841 


11.02 


0.89 


0.96 ± 0.04 


10.58 


1.26 


1.50 ± 0.05 


13.7 ± 


0.8 


190.9 ± 


2.2 


0.19 


NGC 2903 (outer) 


10.45 


1.23 


2.01 ± 0.23 


9.46 


1.78 


1.37 ± 0.85 


18.7 ± 


1.4 


119.3 ± 


1.4 


0.25 


NGC 2976 


8.88 


0.23 


0.42 ±0.11 








<0.1 




>500 




1.65 


NGC 3031 


10.92 


0.96 


1.20 ± 0.65 


9.47 


0.23 


0.23 ± 0.27 


9.3 ± 


60.4 


77.1 ± 


96.5 


3.31 


NGC 3198(1 comp) 


10.27 


0.60 


0.75 ±0.14 








10.1 ± 


1.8 


107.6 ± 


2.8 


1.30 


NGC 3198 (2 comp) 


10.29 


0.55 


0.69 ± 0.03 




0.00» 


0.00=" 


11.2 ± 


0.43 


104.0 ± 


0.7 


0.37 


IC 2574 




0.00" 


0.00" 








<0.1 




>500 




1.66 


NGC 3521 


10.87 


0.43 


0.59 ± 0.21 








14.0 ± 


12.6 


122.5 ± 


20.4 


5.48 


NGC 3621 


10.18 


0.47 


0.80 ± 0.03 








6.5 ± 


0.4 


128.0 ± 


3.3 


0.52 


NGC 4736 


10.09 


0.42 


0.67 ±0.18 


9.51 


0.15 


0.28 ± 0.08 


72.3 ± 


54.9 


43.4 ± 


4.6 


1.44 


DDO 154 


8.00 


1.21 


3.79 ± 2.84 








3.4 ± 


1.3 


68.0 ± 


16.8 


0.81 


NGC 5055 




O.OO" 


0.00" 


9.60 


0.21 


0.16 ±0.01 


33.3 ± 


0.6 


121.6 ± 


0.5 


0.57 


NGC 6946 


10.68 


0.52 


0.82 ± 0.06 


9.35 


0.59 


0.59 ± 0.05 


3.2 ± 


2.6 


281.4 ± 


183.2 


1.01 


NGC 7331 (grad) 


10.90 


0.34 


0.35 ± 0.02 


10.13 


0.75 


0.69 ± 0.07 


9.3 ± 


1.2 


171.2 ± 


8.1 


0.27 


NGC 7331 (const) 


11.08 


0.64 


0.64 ± 0.03 


10.10 


0.90 


0.90 ± 0.09 


3.3 ± 


1.1 


257.3 ± 


56.9 


0.24 


NGC 7793 


8.22 


0.02 


0.06 ± 0.96 








15.3 ± 


9.6 


89.5 ± 


20.7 


3.53 


NGC 7793 (rising) 


8.96 


0.10 


0.33 ± 0.34 








<0.1 




>500 




1.30 



Note. — (2): Logarithm of the predicted stellar mass of the disk (Mq). (3): Average stellar mass-to-Hght ratio in the 3.6 fim band (A/q /Lq). (4): 
Scaling factor with respect to the predicted T^ '^^ value listed in Table 3. The factor with respect to the Kroupa T^ '"' value listed in Table 4 can be 
obtained by multiplying with a factor of 1 .41. (5): Logarithm of the predicted stellar mass of the bulge (Mq). (6): Average stellar mass-to-light ratio 
in the 3.6 fim band (Mq /Lq ). (7): Scaling factor with respect to the predicted ""^ value listed in Table 3. The factor with respect to the Kroupa 
value listed in Table 4 can be obtained by multiplying with a factor of 1.41. (8) NFW parameter c and associated uncertainty (kpc). (9) NFW 
parameter V200 and associated uncertainty (km s~'). (10) Reduced xl- 
^ The stellar mass-to-light ratio T^^ was fixed to T^"^ = 0. 
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Fig.|57]all derived rotation curves on the same physical scale, 
with the position of the origin of the rotation curve offset us- 
ing the absolute magnitude Mb of the galaxy. The well-known 
change i n properties o f the rotation curves with luminosity 
(see, e.g.. iBroeilsll 1 992h is readily apparent. Luminous galax- 
ies have rotation curves that rise steeply, followed by a decline 
and an asymptotic approach to the flat outer part of the curve; 
low-luminosity galaxies show a more gradual increase, never 
quite reaching the flat part of the curve over the extent of their 
H I disks. 

This is shown more clearly in the right hand panel in 
Fig. |57] where we plot the logarithmic rotation curves, log(y) 
vs log{R). Once again the starting point of the rotation curve 
is determined by the absolute luminosity. The change in slope 
as a function of luminosity can readily be made out. 

Another well-known relation for disk galaxies, is that low- 
luminosity galaxies (late-type disk and dwarf galaxies) are 
more dominated by dark matter than luminous galax ies (early- 
type disk galaxies) (e.g.. [Roberts & Haynesl[l994h . The flat 
rotation curves of spiral galaxies are additionally an indica- 
tion that dark matter should become more and more dominant 
towards larger radii. We can explicitly evaluate these relations 
by deriving the ratio of the mass of the baryons and the mass 
of the dark matter halo at each radius as 



Mb. 



aryons 



y2 +T,y2 



obs 



(11) 



where Mbaryons/Mtot is an indication for the degree of dark 
matter dominance, Vobs is the observed rotation curve, Vgas the 
rotation curve of the gas component as described in Sect. 5.1, 
and ^/TlVi, is the rotation velocity of the stellar component 
(taken to be the quadratic sum of multiple stellar components 
if present). 

Note that for some of the galaxies discussed in Sect. 6, we 
derived multiple mass models with, e.g., single or multiple 
stellar components. In the rest of the paper, we will, so as not 
to needlessly confuse the discussion, restrict ourselves to one 
model per galaxy. For NGC 2403 and NGC 3198 we use the 
two-component model, for NGC 7793 the model that uses the 
entire rotation curve, and for NGC 7331 the model without 
the color gradient. These particular choices do not affect any 
of our conclusions. 

The top panel of Fig. |58] shows the radial dependence of 
Mbaryons/Mtot for our Sample galaxies, assuming values for 
T^-* as predicted from the photometry. These distributions are 
for all practical purposes independent of the choice of IMF. 
The largest variation in Mbaryons/Mtot is about 0.1 dex (a fac- 
tor 1.3) in the very inner parts of the luminous galaxies, and 
quickly becomes totally negligible towards lower luminosities 
and large radii. 

There are a few things to note in Fig. |58] Firstly, galax- 
ies with the highest absolute luminosity are in general the 
least dominated by dark matter. This is not only true glob- 
ally, but also locally: even at the outermost radii, the baryons 
in these luminous galaxies are more dominant than in the 
least luminous galaxies. Secondly, while the most lumi- 
nous galaxies show a clear gradient in Mbaryons/Mtot with ra- 
dius (with the outer radii becoming progressively more dark- 
matter-dominated), this trend is less pronounced in the low- 
luminosity galaxies. The latter indeed are dominated by dark 
matter at all radii. These results are robust against any par- 
ticular choice of dark matter model, but only depend on the 
values of Tj^. 



It is also interesting to explore the distribution of 
Mbaryons/Mtot for the casc of TJ-^ as a free parameter (i.e. as 
determined by the dynamics, rather than the photometry). In 
these cases the results are not entirely independent of the halo 
model assumed, as the choice of model plays a role in deter- 
mining T^^ during the fitting. The middle and bottom pan- 
els of Fig. |58]show Mbaiy-ons/Mtot for the best-fitting ISO and 
NFW T^ '' values, respectively. The trends discussed above 
are still present in the fixed T^ *" ISO case, but they become 
confused for the NFW model. The latter is most likely a re- 
flection of the NFW model being an inappropriate fitting func- 
tion for at least some of the galaxies. 

7.2. Stellar mass-to-light ratio trends 

An overview of the derived T^-* values of the main disk 
components of our sample galaxies is given in Fig.|59l Here 
we plot the photometric (fixed) and dynamical (free) values 
for TJ-^ against the colors and luminosities of the galaxies. 
The left-hand panels show the fixed and free values of 
the disks for both the ISO and NFW model plotted against 
(J-K) color We also distinghuish between our two choices 
for the IMF. Also shown is the predicted relation between Tj^ 
and color (Eqs. 4 and 5), again for both the diet-Salpeter and 
the Kroupa IMFs. By definition, the photometric Tj '' values 
follow their respective relations. The small scatter is caused 
by the color gradients present in these galaxies which make 
the photometric T^-* differ slightly from those with a constant 
T^-* as a function of radius. 

In general, we see that the diet-Salpeter curve defines an 
approximate upper limit to the distribution of the majority 
of the best fit values (apart from a few obviously discrepant 
cases). Accepting the best fit (free) T^'^ values at face value, 
this would suggest that a diet-Salpeter T^^ analysis overesti- 
mates the disk masses slightly. Indeed, the Kroupa fixed T^^ 
values seem to be a better match to the free Tj^ values. Alter- 
native explanations could be effects of star formation history, 
or unexpectedly large contamination by PAHs, AGBs or hot 
dust in the 3.6 /im maps, though it is likely that the IMF plays 
the dominant role. Future rigorous population synthesis mod- 
eling should shed light on some of these issues. The NFW 
T^-^ values tend to be lower than the ISO values, which is to 
be expected due to the intrinsically steeper mass distribution 
of the NFW halo. 

We compare the distributions of the T^-^ values in Fig.l60l 
The diet-Salpeter fixed Tj'^ values are offset towards some- 
what higher values compared to the ISO free T distribution. 
The Kroupa fixed Tj-^ values are a much better match, as al- 
ready suggested by Fig. |59] The distribution for the NFW 
model is less well-defined, and has a larger number of galax- 
ies at very low T^-^ values. These are the cases where T^^ 
had to be set to zero manually, due to the preferred negative 
T^ *" values. The smaller spread of the T^^ values for the ISO 
case, and the somewhat better agreement with the predicted 
values, especially for the Kroupa IMF results, could be taken 
as tentative supporting evidence that the ISO model is supe- 
rior in describing the properties of real galaxies. 

In general, it is remarkable that despite our empirical and, 
from the point of view of a population synthesis expert, un- 
doubtedly rudimentary modeling of the stellar mass-to-light 
ratios, the best-fit dynamical T;^ values lie so close to the 
photometric T^^ values. The photometric Tj^ values "don't 
know" anything about dynamics, and the best-fitting dynami- 
cal values, whereas the best-fitting T^^ ratios have no a priori 
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Fig. 3 1 . — ISO and NFW rotation curve fits for NGC 2903. Only radii R > 3.3 kpc are considered liere. Lines and symbols as in Fig. 1261 The horizontal line 
in the top-left panel represents the degenerate halo rotation curve; see text for a description. 



reason to stay close to the predicted photometric stellar pop- 
ulation values. It is therefore gratifying to see that the calcu- 
lated TI'^ agree with the fitted T^^ to such a good degree. 

Plotting TJ-^ against luminosity (right-hand panels in 
Fig. l59l ) we see that the brightest galaxies have the highest 
T^-^ ratios. These large values are not found for the lower 
luminosity galaxies. 

The one galaxy that shows a different behavior is DDO 154. 
Its predicted population-synthesis Tj^ value is low, but the 
T^-^ values derived from the mass-modeling do not agree and 
are much higher This discrepancy can neither be due to the 
data quality (THINGS observed galaxies at identical resolu- 
tion and sensitivity), nor due to large non-circular motions, 
as DDO 154 has a well- behaved, "quiescent" velocity field 
jTrachternach et al.ll2008l) . 

One can interpret this result at face value, and assume that 



DDO 154 does have a massive disk. This would, however, 
put it out of odds with the rest of the sample, and imply a sub- 
stantial change in the nature of galaxy disks towards low lumi- 
nosities. An alternative explanation could be that for galaxies 
like DDO 154 the disk is not massive enough to have a notice- 
able effect on the total rotation curve. There would then be no 
small deviations induced by, e.g., spiral arms, and the dynam- 
ics are everywhere dominated by the dark matter. The fitting 
program is then unable to derive a sensible value for T^ ''. 
Obviously, a larger sample of high-resolution observations of 
galaxies fainter than Mb ~ -17 should then show the same 
trend. A similar, but slightly more prosaic, explanation, is re- 
lated to the fact that DDO 154 has the lowest surface bright- 
ness disk in 3.6 i-im in our sample. The surface brightness 
profile is only measured out to ^ 100" 2 kpc), while the ro- 
tation curve extends out to ~ 400" (~ 8 kpc). This means that 
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Fig. 32. — Surface brightness and T^^ profiles for NGC 2976. Lines and symbols as in Fig.l24l 



most of the disk rotation curve is based on a smooth exponen- 
tial extrapolation of the observed surface brightness profile. 
If the underlying (undetected) disk does indeed induce small 
wiggles in the observed rotation curve, then these will not be 
matched by the disk rotation curve (as derived from the expo- 
nential extrapolation), and this will throw the fitting program 
off track. Tests of replacing the observed disk rotation curve 
of a few of the other THINGS galaxies with their smooth ex- 
ponential disk equivalents, indicate that this does indeed af- 
fect the derived Tj^ values, but the effects are not unambigu- 
ous enough to conclusively show that this is the single cause 
of what is happening in the DDO 154 mass modeling. Deeper 
observations of the stellar disk should help pinpoint the exact 
cause of the discrepancy. 

Finally, in Fig. [61] we show the dynamical importance of 
the stellar disk as a function of luminosity. We plot the ratio 



of the maximum rotation velocity of the disk (assuming the 
photometric T^ *" ratios and using both IMF assumptions) and 
the total rotation velocity at the radius where the maximum 
disk velocity occurs (2.2/2 for an exponential disk). It is clear 
that disks become dynamically less and less important with 
decreasing luminosity. Note that a similar plot of the impor- 
tance of the disk versus luminosity for the case of free T^^ 
(also shown in Fig.l6TI) does not show such a trend as clearly, 
but here the result may b e diluted by the e ffects just described. 

In a previous analysis. iBottemal (1 1997h used measurements 
of stellar velocity dispersions in a number of early-type disk 
galaxies to estimate the dynamical contribution of the disk. 
He found for the magnitude range he investigated (-22 < 
Mg < -19) that the ratio of maximum rotation velocity of the 
stellar disk and total rotation velocity at the corresponding ra- 
dius (generally 2.2h) is VdtkV^obs'' = 0.63 ±0.10. These re- 
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Fig. 33. — ISO and NFW rotation curve fits for NGC 2976. Lines and symbols as in Fig. 125 1 



suits are also indicated in Fig. 1611 

Over the same range in magnitude we find that V^^^^^ /V^^f' = 
0.81 ± 0.19 for the fixed diet-Salpeter T^*" case and 
^diskV^obs'' = 0-68 ± 0.16 for the Kroupa IMF fits. The latter 
is in reasonable agreement with the Bottema (1997) values. 
The free T^ '' results yield VdikV^obs'' = 0-73 ± 0.17. 

In order to gauge whether our data show a preference for 
one particular IMF, we compare in Fig.|62]the goodness-of-fit 
values for the fixed T^^ models (for both the diet-Salpeter and 
Kroupa IMFs) and the ISO and NFW halos as a function of lu- 
minosity. Several conclusions can be drawn from this compar- 
ison. Firstly, the choice of IMF does not influence the good- 
ness of fit for low-luminosity galaxies. This is of course to be 
expected as the dynamics of these galaxies are dominated by 
the dark matter component. Secondly, for the intermediate lu- 
minosity galaxies the ISO models prefer a diet-Salpeter IMF, 



while the NFW models prefer a Kroupa IMF. This is also eas- 
ily understood, as the NFW halo has a steeper mass-density 
profile (a steeper inner rotation curve) which leaves less room 
for the luminous component and results in lower optimal Tj^ 
values. Thirdly, for the high-luminosity galaxies, the Kroupa 
IMF is much preferred over the diet-Salpeter fits. This reflects 
the fact that the diet-Salpeter fits result in disks that are close 
to the maximum allowed by the dynamics. A small change 
in the Tj-^ value can then result in a large difference in the 
quality of the fit. The Kroupa IMF with its lower T^^ values 
is more stable against these effects, as it still allows a dark 
matter component in the inner parts which can compensate 
for any small Tj^ uncertainties. 

The photometrically determined, fixed T^^ values there- 
fore do not allow us to choose a particular IMF, without mak- 
ing additional assumptions on the distribution of the dark mat- 
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Fig. 34. — Surface brightness and T^^ profiles for NGC 3031. Lines and symbols as in Fig. 1241 



ter. It would thus be interesting to independently confirm the 
values of the disks with measurements of the stellar veloc- 
ity dispersions within the THINGS disks, and thus calibrate 
our measurements. 

7.3. Halo rotation curves 

In Sect. l7.1l we used the derived T^^ values in combination 
with our rotation curves to derive the ratio of baryonic to dark 
matter as a function of radius. In a similar way we can derive 
the rotation curves of the dark matter halos: 

where again we are, strictly speaking, dealing with the cir- 
cular rotation velocities of test particles in the plane of the 
galaxy. The velocities of the gas and stars have been cal- 
culated using the three-dimensional mass distributions as de- 
scribed in Sect. 5.1 and 5.2. We again assume spherical halos. 



iMcGaugh et al.l (l2007h show that these empirical halo rota- 
tion curves, for a realistic range of values, yield halo dark 
matter densities about ~ 50 percent lower than expected from 
ACDM simulations. One interesting aspect of their conclu- 
sions is that they are based on the outer (R> \ kpc) parts of 
the rotation curves, i.e., those regions where the usual con- 
cerns about resolution or slope of the curve are not relevant. 

Here we repeat the analysis of iMcGaugh et al.l ( l2007h for 
our current sample. We derive the halo rotation curves for our 
galaxies assuming the photometric (fixed) T^'^ values for both 
IMF assumptions. The resulting diet-Salpeter IMF logarith- 
mic halo rotation curves are shown in Fig.|63] where we have 
divided the galaxies into three bins according to their lumi- 
nosity. The Kroupa IMF halo rotation curves are very similar 
and we will not show them here, though we do discuss their 
implications below. 
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ISO and NFW rotation curve fits for NGC 3031. Lines and symbols as in Fig. 1261 



Fig. 35.— 

For each luminosity bin we show the original data, as well 
as the observed average halo rotat ion velocity, binned in steps 
of 0.25 dex in radius. Following iMcGaugh et alJ (l2007h . we 
only consider radii R > I kpc, so as not to confuse the is- 
sue with potential systematic effects in the determination of 
the very inner slope. Over the radial range observed, the av- 
erage velocities behave very nearly like a power-law. Also 
shown in Fig. |63] are two NFW model rotation curves with 
V2(K) velocity values chosen so they encompass the observed 
data set. The lowest velocity model has V200 = 50 kms"', 
and is representative of the dwarf galaxies in our sample. 
On the high-luminosity end of the range we show a model 
with V200 = 300 kms"'. Galaxies with intermediate luminosi- 
ties or masses will have curves lying between these two ex- 
tremes. The corresponding c-parameters of these curves have 
been chosen to be consistent with current ACDM cosmologi- 



cal models (cf. Sect. 5.3.1 and 7.4.2). 

In Fig.|63]we also show the fit to the halo rotation ve l ocities 
derived for the ~ 60 galaxies in the iMcGaugh et alJ (l2007h 
sample, assuming their photometrically derived ratios. 
They find that over the radial range observed, the halo rotation 
velocity can best be describe by log V = 0.50 log 1 . 49. To 
compare this with CDM predictions, IMcGaugh et al.l (l2007h 
subjected a large number of NFW halo models to observa- 
tional selection effects and found that the average halo rota- 
tion velocity in a ACDM universe is best described as log V = 
0.48 log R+1.66 (also shown in Fig.l63]l. From this they con- 
clude that the empirically determined densities of dark matter 
halos are ~ 50 percent lower than those predicted by CDM 
simulations. 

We first consider the lowest luminosity bin with Mb > 
-18.5. Over the observed radial range the average diet- 
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Fig. 36. — Surface brightness and T^^ profiles for NGC 3198. Lines and symbols as in Fig. 1241 



Salpeter halo velocities can be described by a power-law 
log V = 0.52 log /?+ 1 .25 (Kroupa: log V = 0.41 log /?+ 1 .35). 
The normalization is thus 0.24 (0.15) dex (a factor 1.7 to 
1.4) lo wer than the empirical one derived bv lMcGaugh et al] 
(120071) and 0.41 (0.31) (a factor 2.5 to 2.0) dex lower than the 
average velocity for CDM halos as derived in that paper. The 
latter value was, however, derived for a much larger range 
in mass than considered in this low-luminosity sub-sample. 
Nevertheless, even if we compare with the value found at 1 
kpc for a V200 = 50 kms"' halo (which should be more typical 
for the galaxies considered here; see Fig. |63] ). we still find a 
velocity that is 0.25 (0.15) dex lower than predicted. Only at 
R ^ 10 kpc do we start finding some agreement between the 
empirical and simulated velocities. 

The situation is similar for the intermediate luminosity 
galaxies with -20.5 < Mb < -18.5. The average diet-Salpeter 



derived halo velocity can be described as log V = 0.50 log R+ 
1 .45 (Kroupa: log V = 0.47 log /?+ 1.50) , virtually indistin- 
guishable from the lMcGaugh et all (l2007h fit. The same con- 
clusions thus hold here: the velocity normalization is about 
0.21 (0.16) dex (a factor 1.6 to 1.4) too low compared to aver- 
age CDM halos, or about 0.25 (0.20) dex (a factor 1.8 to 1.6) 
lower than that of a V200 =150 kms~' halo (which has a value 
logV ~ 1 .7 at 7J = 1 kpc). This difference in normalization is 
illustrated by the fact that in Fig. |63]the inner three observed 
average data points overlap with those of the V200 = 50 kms"' 
halo, even though the average maximum rotation velocity of 
the galaxies in this sub-sample is ~ 130 kms"'. 

At the high-luminosity end of the sample, Mg < -20.5, 
all observed average data points are now within the range 
described by NFW halos, but there is still a normalization 
discrepancy for the diet-Salpeter case. The data extend fur- 



50 



de Blok et al. 



150 



100 

E 



> 



50 







150 



100 

E 



> 



50 - 







— I — I — I — I — I — r 

ISO, fixed 




'> \-''\ 1 












'> \-''\ 1 










1 1 1 1 


1 1 1 1 


1 1 1 1 








1 1 1 1 


1 1 1 1 


1 1 1 1 



NFW. fixed 




J I I I I I I I I I I I I L_ 



ISO, free 



NFW, free 




1- I-' I I I I I I I I I I I I I I I I I I i 



10 20 30 

R (kpc) 



40 10 20 30 

R (kpc) 



40 



Fig. 37. — ISO and NFW rotation curve fits for the multiple-disJc model of NGC 3198. Lines and symbols as in Fig. 1261 



ther out in radius and are best described by a power-law fit 
log V = 0.32 log /?+ 1 .69 (Kroupa: log V = 0.27 log /?+ 1 .83). 
The slope is thus more shallow than in the other luminosity 
bins, perhaps indicating we are starting to see the velocity 
turn-over or flattening expected in the ISO and NFW mod- 
els. If we look at halos with V^oo values between 200 and 
300 kms"', we find intercept values at = 1 kpc between 
logV = 1.75 and 1.80. This discrepancy is therefore less pro- 
nounced than for the lower luminosity bins, but still amounts 
to ~ 0.1 dex, or a factor 1.3 difference, for the diet-Salpeter 
case. For the Kroupa case, the intercepts agree, but note that 
the slope of the average halo rotation curve is much shallower 
than the NFW model curves. The real rotation curves flatten 
much faster than NFW models with comparable mass. 

For completeness we note that the halo velocities of the en- 
tire THINGS rotation curve sample (assuming a diet-Salpeter 



IMF) are best descri bed as log V = . 46 log R + 1.49, in good 
agreement with the iMcGaugh et alj ( |2007|) result. For the 
Kroupa case we find log V = 0.45 logR+ 1.58. The best-fit 
dynamical, free T^^ values (as derived for the ISO model) 
give a fit log V = 0.50 log R+l .50, indistinguishable from the 
fixed Tl^ diet-Salpeter IMF result. 

In summary, the empirically derived rotation velocities 
associated with dark matter halos of galaxies fainter than 
Mb ~ -19 are about 0.2 dex or ~ 1.6 times lower than pre- 
dicted by CDM simulations, independent of IMF, with the 
largest discrepancies found in the lowest luminosity galaxies. 
The agreement with simulated CDM halos improves towards 
higher luminosities, but the degree of improvement strongly 
depends on the assumed IMF. 

If we translate the 0.2 dex in difference in velocity to den- 
sity differences, we find an average density that is ^ 40 to 
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Fig. 38. — ISO and NFW rotation curve fits for the single-disk model of NGC 3198. Lines and symbols as in Fig. 1251 



~ 60 percent of the predicted values. These numbers do de- 
pend on the precise density distribution, but they suggest that 
dark matter halos of late-type galaxies have densities that are 
about half of those predicted by CDM. This may also hold 
true of the luminous galaxies in our sample, but here the re- 
sults depend on the assumed IMF. 

7.4. Mass models and model parameters 
1 A.\. Treatment of uncertainties 

In the preceding discussion we have analyzed the impor- 
tance of dark matter in the THINGS galaxies without any 
strong assumptions on its actual distribution. However, ro- 
tation curves are most often used to test particular dark matter 
models, and, as described in Sect. 1, the debate in recent years 
has focused on the apparent inability of cosmological CDM 
simulations to describe the distribution of dark matter on the 



scale of galaxies. With the high resolution and excellent sen- 
sitivity of the THINGS velocity fields and rotation curves, we 
are now in a position to revisit some of the outstanding ques- 
tions, and compare the applicability of both the observation- 
ally motivated ISO model and the cosmologically motivated 
NFW model (cf. Sect. 5.3). 

Comparison of the fits derived here can tell us which of the 
dark matter models, if any, is preferred. Note that with the 
sampling of our velocity fields of two points per beam, the 
individual rotation curve data points are not entirely indepen- 
dent. This obviously affects the absolute values, but with 
the non-Gaussian distribution of the uncertainties in the rota- 
tion curve data points, we will not dwell on an interpretation 
of the absolute values. More importantly, our conserva- 
tive estimates of the uncertainties, taking into account the dis- 
persion of velocity values found along tilted rings, as well 
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Fig. 39. — Surface brightness and T^'^ profiles for NGC 3521. Lines and symbols as in Fig.l24l 



as the differences between approaching and receding sides 
(Sect. [l!6] l will also affect the absolute x" values. More rele- 
vant for our discussion are therefore the differences between 
the reduced x" values, which tell us the relative merits of any 
particular model. In this respect, the -statistic should not be 
confused with the x^-test. The statistic simply quantifies the 
sum of the absolute differences between the observed veloc- 
ities and the best-fitting model velocities, expressed in terms 
of the measurement uncertainties and normalized using the 
effective number of data points. This reduced is here only 
used as a proxy for the average excess velocity per data point. 

In determining the halo parameters and the correspond- 
ing uncertainties the ROTMAS software uses a simple least- 
squares procedure. The uncertainties in the rotation velocity 
should therefore be reflected in the uncertainties in the halo 
parameters, and our definition of the rotation velocity uncer- 



tainties could thus potentially have an impact on the values 
we derive. We evaluated the effect of adopting a different def- 
inition for the uncertainties (as discussed in Sect. 13.6b on the 
fit parameters, and found the resulting differences to be small, 
and certainly negligible when compared to the differences in- 
troduced by small variations in T\ . We illustrate this for two 
typical galaxies in Table|7] where we Ust the fitted parameters 
for the ISO and NFW halos, assuming fixed T^ '' values, and 
adopting respectively, f i) our original definition for the uncer- 
tainties, {ii) the dispersion of velocities along the rings, and 
{ Hi) the formal fit error in the rotation velocity. 

The extreme values for as found when adopting the 
formal fit error from ROTCUR, immediately show that these 
severely underestimate the true uncertainties, as already dis- 
cussed in Sect. 13. 61 It is clear that for the other two definitions 
of the errors small differences in the resulting halo parame- 
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Fig. 40. — ISO and NFW rotation curve fits for NGC 3521. Lines and symbols as in Fig.l25] 



ters exist, but it is also clear that different (reasonable) defi- 
nitions of the uncertainties do not dominate the error budget: 
the largest uncertainty is introduced by the stellar mass-to- 
light-ratio T^, (cf. the results in Table |7] with those in Tables 
[3]|6]). In summary, our definition and treatment of the uncer- 
tainties in the rotation velocities does not dominate the final 
uncertainties in the halo parameters. 

7.4.2. Halo parameters 

We now consider the decompositions with T^^ as a free 
parameter (Tables |5] and |6l). The corresponding values are 
compared in Fig.|Ml The left panel compares the values of x^- 
for both models, plotted against absolute luminosity. Three 
galaxies (xj. > 2) are obviously not well described with either 
model. Concentrating on the other galaxies, we see a sys- 
tematic difference between high- and low-luminosity galax- 



ies. The luminous galaxies show comparable values for 
either model; the low-luminosity galaxies with Mb > -19 are 
clearly better described using the ISO model. This is shown 
more clearly in the right hand panel of Fig. |64l where we 
show the difference between the xl values for both models 
as a function of luminosity. 

NFW and ISO models thus fit equally well for high lumi- 
nosity galaxies, with reasonable T^^ values. The extent to 
which this is determined by subtle trade-offs between, e.g., 
dynamical contributions of bulge, disk and halo still needs 
to be determined, but such an analysis is likely limited by 
the accuracy with which stellar mass-to-light ratios can be 
determined from photometric data and models. Alternative 
methods (via, e.g., dispersion measurements) will be needed 
to constrain the importance of the disk fur ther. 

Regarding the low-luminosity galaxies, iTrachternach et all 
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Fig .41. — Surface brightness and T^'^ profiles for NGC 362 1 . Lines and symbols as in Fig.|24] 



(Hogs') show, & om harmonic decompositions of the velocity 
fields, that non-circular motions are lowest in these galax- 
ies. The increasing incompatibility of the NFW model model 
towards lower luminosities, is therefore difficult to explain 
purely in terms of these non-circular motions. Note also that 
the current sample has been observed at identical and homo- 
geneous sensitivity and resolution and has been analyzed in a 
consistent manner The differences found here can therefore 
not be attributed to different data qualities or analysis meth- 
ods. 

Lastly, we investigate the model parameters. To start with 
the NFW model, as described before, the fit parameters c and 
V200 are related and their valu es are determined by t he as- 
sumed cosmology (Sect 5.3.1). iMcGaugh et alJ (l2007h show 
the relation between c and V200 for the "vanilla ACDM" model 
presented in lTegmark et al.l (l2004l) . and the 3-year WMAP re- 



sults (S^ergeretal] (l2007h . Both show a fairly flat distribution 
of c as a function of V200, but with a different offset (due to 
a different power spectrum normalization). Both distributions 
are shown in Fig. l65]along with the la spread expected from 
simulations (Bulloc k et al.ll2001h . It is clear that for rotation 
curves which are consistent with CDM, we expect to find con- 
centration values between c ^ 4- and c ~ 10, fairly indepen- 
dent of Vaoo- In Fig. |65] we also plot the values for c and V200 
derived from our rotation curve fits, both for photometrical 
(diet-Salpeter and Kroupa IMF) and dynamical T^ '' values. 
Also indicated are the galaxies for which no sensible fitting 
parameters could be derived. It is immediately clear that the 
distribution of c values we find shows much more scatter than 
the CDM relations. At first glance one could conclude that we 
are finding c-values that are on average higher than predicted, 
but it is not clear to what extent our distribution is affected by 
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Fig. 42. — ISO and NFW rotation curve fits for NGC 3621. Lines and symbols as in Fig. 125 1 



the absence of galaxies for which sensible NFW parameters 
could not be derived. 

As the ISO halo has no basis in current cosmology, 
there are no a priori expectations for its model param- 
eters. One can however still derive them and check 
for poss ible trends with other f undam ental galaxy parame- 
ters, as iKormendv & FreemanI (120041) have done. Fig. |66] 
shows the distribution of the central density po against the 
core radius Rc. The distr ibution is consistent with the 
iKormendy & FreemanI (l2004l) relation (their Eq. 20), but with 
slightly larger scatter 

Similar relations with the absolute luminosity (shown in 
Fig. Wh are less well-defined. The most noteworthy aspect 
is that most galaxies have a core radius Rc> I kpc, indepen- 
dent of the adopted value for Tj^, showing that the effects of 
the "soft" core are already clearly noticeable well outside the 



central regions. In other words, while scales smaller than ~ 1 
kpc are needed for a clear distinction between cusp and core, 
"symptoms" of the core should already be visible at scales 
larger than that. 

8. SUMMARY 

We present a uniquely high-resolution analysis of the kine- 
matics of 19 galaxies. The galaxies formpart of THINGS, and 
have therefore been observed and analyzed in a homogeneous 
and consistent manner The high spatial and velocity resolu- 
tion, as well as the exquisite sensitivity enable us to derive the 
tilted-ring parameters based solely on the H 1 data. There is 
therefore no dependence on, e.g., optical inclinations or axis 
ratios, which generally carry their own systematic effects with 
them. We summarize our results as follows: 

• We determine the 3.6 /im stellar mass-to-Ught ratio T^^ 
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Fig. 43. — Surface brightness and T^^ profiles for NGC 4736. Lines and symbols as in Fig. 1241 



in two fully independent w ays: one is bas ed on pop- 
ulation synthesis modeling jOh et al.ll2008l) . the other 
is a purely 'dynamical' Tj-*, determined solely by the 
best fit T^ '' value in the rotation curve mass model. 
For the population synthesis Tj-^ values we consider 
two assumptions for t he IMF: the "diet"-Salpeter IMF 
(iBell & de Jong 1200 Ih. design ed to fit maximum disk 
models, and the lKroupal (1200 lb IMF, which yields disks 
that are ^ 40 percent lighter We find good agreement 
between the photometric and dynamical Tj^ values, 
especially for the Kroupa IMF. The close match be- 
tween the dynamical and population synthesis results 
indicates the stellar disks of the galaxies in the THINGS 
sample do not contain large amounts of dark matter as- 
sociated with the disks themselves. The dark matter 
content of the THINGS galaxies is associated with the 



halo, not the disk. Our analysis also obviates the need 
for so-called minimum or maximum disk approaches. 

• We do not find any steeply declining rotation curves out 
to the last measured point. Specifically, previous ob- 
servations seemed to indicate declining rotation curves 
for NGC 3521, NGC 7793, NGC 2366 and DDO 154. 
With the increased resolution of the data we can now 
better trace the inclination and PA trends. Any gentle 
decreases in the velocity that we still find can now be 
attributed to uncertainties in the rotation velocity or in- 
clination. There is therefore no evidence from THINGS 
that we have probed the "edge" of a dark matter halo. 

• We fit NFW and ISO models to our rotation curves us- 
ing our best estimates for T^^. We find that for galax- 
ies Mb < -19 both models statistically fit equally well. 
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Fig. 44. — ISO and NFW rotation curve fits for NGC 4736. Lines and symbols as in Fig.l26l 



Note, however, that for the 11 galaxies in this high- 
luminosity bin, the fixed T^^ diet-Salpeter NFW fits re- 
sult in 4 galaxies with unphysical (severely inconsistent 
with CDM) halo parameters. These unphysical cases 
disappear when a Kroupa IMF is assumed. However, 
the resulting halo parameters still show a scatter much 
larger than the simulation predict. 

• For galaxies with Mb > -19 the core-dominated ISO 
model fits significantly better than the NFW model 
(and the fixed NFW fits yield 4 out of 5 galax- 
ies with unphysical halo parameters independent of 
IMF assumption). The success of the ISO halo for 
the se galaxies is remarkabl e, given the conclusion 
by iTrachternach et all (|2008|) that the THINGS low- 
luminosity galaxies have the smallest non-circular mo- 



tions: the observed dynamics most likely reflect that of 
the dark matter halo. 



• We find that the low-luminosity galaxies in our sample 
are dark matter-dominated throughout, based on both 
the population synthesis Tj-^ values, as well as the 'dy- 
namical' T^-^ results. As noted above, the agreement 
between these two sets of values indicates the dark mat- 
ter is not associated with the stellar disk. The dark mat- 
ter rotation curves of these galaxies are about 0.2 dex 
(a factor 1.6) lower than predicted by CDM over their 
entire disks. This implies that the density of dark mat- 
ter in real late-type galaxies is ^ 50 percent of what is 
predicted by CDM simulations. 
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Fig. 45. — Surface brightness and T^'^ profiles for DDO 154. Lines and symbols as in Fig. [24] 
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Fig. 46. — ISO and NFW rotation curve fits for DDO 154. Lines and symbols as in Fig.lls] 



TABLE 7 

Uncertainties in model parameters 



ISO halo NFW halo 



Name 


type 


A 


Rc 


PO 




c 


1^200 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 2403 


total 


4.64 ±1.09 


2.14±0.05 


77.9 ±3.3 


0.79 


9.8 ±0.2 


110.2±1.0 


0.56 




disp 


4.22 ±1.21 


2.16±0.05 


77.1 ±3.3 


0.94 


9.9 ±0.2 


109.5 ±0.9 


0.70 




fit 


0.21±0.15 


2.07 ±0.05 


80.2 ±3.0 


399.24 


10.5 ±0.2 


105.2± 1.2 


403.93 


DDO 154 


total 


2.40 ±0.91 


1.34±0.06 


27.6 ± 1.6 


0.44 


4.4 ±0.4 


58.7 ±4.3 


0.82 




disp 


2.25 ±0.90 


1.31 ±0.05 


28.4± 1.7 


0.51 


4.4 ±0.4 


59.0±4.2 


0.87 




fit 


0.13 ±0.07 


1.44 ±0.06 


24.7 ± 1.5 


141.15 


5.6±0.3 


49.5 ±2.0 


168.37 



Note. — (1) Name of galaxy; (2) definition of uncertainties; "total": definition adopted in this paper; "disp": dispersion in 
velocities along tilted rings; "fit": formal fit eiTor in rotation velocity; (3) Average value and dispersion of uncertainties; (4)-{9) halo 
fit parameters, definitions as in Tables f3T6| 
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Fig. 47. — Surface brightness and T^'* profiles for NGC 5055. Lines and symbols as in Fig.|24] 
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Fig. 48. — ISO and NFW rotation curve fits for NGC 5055. Lines and symbols as in Fig.l26l 



62 



de Blok et al. 




Fig. 49. — Surface brightness and T^'* profiles for NGC 6946. Lines and symbols as in Fig.|24] 
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Fig. 50. — ISO and NFW rotation curve fits for NGC 6946. Lines and symbols as in Fig.l26l 
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Fig. 52. — ISO and NFW rotation curve fits for the color-gradient model of NGC 7331. Lines and symbols as in Fig. 1261 
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Fig. 53. — ISO and NFW rotation curve fits for the constant * model of NGC 7331. Lines and symbols as in Fig. 1261 
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Fig. 55. — ISO and NFW rotation curve fits for NGC 7793 using the entire observed extent. Lines and symbols as in Fig. 1251 
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Fig. 56. — ISO and NFW rotation curve fits for NGC 7793 using only radii R < 4.7 kpc. Lines and symbols as in Fig. 1251 
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Fig. 57. — All THINGS rotation curves plotted in linear units in the left panel and in logarithmic units in the right panel. The origin of the rotation curves has 
been shifted according to their absolute luminosity as indicated on the horizontal axis. The bar in the respective panels indicates the radial scale. 
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Fig. 58. — The radial distribution of Mbaryons/Mot in the THINGS galaxies for several assumptions on Tj *. The top panel shows the distribution assuming 
the fixed, predicted diet-Salpeter Tj * values. The center and bottom panels show the distributions assuming the best-fitting Tj * values determined by using the 
ISO and NFW models. Light-gray i^ill curves represent galaxies brighter than Mb = -20.5. Dashed dark-gray curves show galaxies with < -20.5 < Mb < -18.5. 
Black curves show galaxies fainter than Mb = -18.5. 
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Fig. 59. — TJ'' plotted against color and luminosity. Results are shown as a function of colour (J — K) in the lefthand column, as a function of absolute 
magnitude Mb in the righthand column. Top row shows results derived for the ISO halo; bottom row those for the NFW halo. Open squares indicate the results 
derived for the free (dynamical) Tj-* fits, the connected grey symbols show the fixed (photometric) Tj^ values. The upper symbols (open grey circles) indicate 
the results for a diet-Salpeter IMF, the lower symbols (filled grey circles) those for a Kroupa IMF. In the lefthand panels the curves shows the theoretical color- 
* relations as derived from Eqs. (4) and (5). The upper curve assumes a diet-Salpeter IMF, the lower curve a Kroupa IMF. In the top row the arrows indicate 
the high free Tj-* value derived for DDO 154. 
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Fig. 60. — Histograms of the various Tj^ values. Top: predicted values: full histogram assvmies a diet-Salpeter IMF, dashed histogram assvmies a Kroupa 
IMF; center: best-fit values assuming the ISO halo model; bottom: best-fit values assuming the NFW model. 
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Fig. 61 . — Ratio of the inferred maximum rotation velocity of the disk and the total rotation velocity occurring at that radius, plotted against absolute luminosity. 
Open squares show the ratios assuming the free (dynamical) Tj^ values. The connected symbols show those derived for the fixed (photometric) T^'^ values. The 
upper symbols (filled black circles) show the results for a diet-Salpeter IMF, the lower (filled gray circles) show those for a Kroupa IMF. The gray box indicates 
the range in luminosity and velocity ratios derived by Bottema ( 1997) on the basis of stellar velocity dispersion measurements. The dotted line represents a 
least-squares fit assuming the predicted T^^ values. The upper line is derived for the diet-Salpeter IMF and decreases as ~ O.lSlogLg. The lower line assumes 
a Kroupa IMF and decreases as ~ 0.201ogLB. These slopes are, however, to a large extent determined by the DDO 154 result at Mb ~ —14. 
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Fig. 62. — Comparison of the quality of the models under various assumptions for IMF and dark matter halo model. The open circles show the difference in 
goodness-of-fit between a diet-Salpeter Tj * model and a Kroupa Tj-* model assuming an NFW halo. The filled circles show the same for the ISO model. 
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Fig. 63. — Halo rotation curves of the THINGS galaxies. The galaxies are binned according to their luminosities, as indicated in the figure. Individual halo 
rotation curve points derived assuming fixed * values are shown as the small lightgray points. The average velocity values binned every 0.25 dex in radius are 
shown as the large black squares. The thick long-dashed dark grey Hne shows the best power-law fit to those points. The thick black dashed curves indicate the 
minimum and maximum rotation velocities one expects to find for NFW halos at each radius. The lower curve has parameters (c,y20o) = (9.5, 50 kms"'). The 
upper curve has (c, V200) = (6.7, 300 km s"' ). All reahstic and CDM-consistent NFW halos a re expected to fa l l betw een these two curves. The lower of the two 
thin dotted line indicates the best power-law fit to the empirical halo rotation curves derived in lMcGaugh et al.l |2007). The upper one shows their best power-law 
fit to a comparable sample of NFW halos. These curves assume a diet-Salpeter IMF; results for the Kroupa IMF are very similar, see text for a more extensive 
description. 
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Fig. 64. — Left: Distribution of as a function of luminosity with Tj * as a free parameter. Filled circles: ISO fits. Open circles NFW fits. Right: Difference 
in ISO and NFW values. Stars: good fits with xi(ISO)< 2 and xl(NFW)< 2. FUled squares: fits with xi(ISO)> 2 or x|(NFW)> 2. 
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Fig. 65. — Left: Distribution of c as a function of 1^200 for the various assumptions on T^'*. Upper panel shows the results assuming the fixed (photometric) 
T^'* values with a diet-Salpeter IMF. The center panel show the same but for a Kroupa IMF. The bottom panels uses the free T^^ values. Arrows in the bottom- 
right comers indicate the fits with c < I. The hatched area bordered by the full curves indicates the predicted c — V200 relation based on the "vanilla" ACDM 
cosmology pres ented in.Tegmark et aU ilOOA) . The counter-hatched ai'ea bordered by the dashed curves shows the relation that can be d erived from the 3-ve ar 
WMAP results jSpergel etall2007l) . The widths of the distributions correspond to the itlo" scatter in c as derived from CDM simulations jBuUock et aU200ll) . 
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Fig. 66. — Left: Distribution of core density p of tlie ISO lialo against core radius Re- Filled circles show the results derived assuming fixed T^^ values, with 
black symbols assuming a diet-Salpeter I MF and grey circles assuming a Kroupa IMF. Open circles show the results derived using free T^^ values. The full 
lines show the relation derived in ,Kormendv & FreemanI 420041) along with the la scatter. 
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Fig. 67. — Left: Distribution of ISO halo core radius Rc (left panel) and core density po (right panel), plotted against absolute luminosity Mb. Filled symbols 
indicate the fixed Tj * results, with black circles representing the diet-Salpeter fits and grey circles the Kroupa fits. The open circles show the free Tj * outcomes. 
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APPENDIX 

DATA AND ROTATION CURVE DERIVATION 

We present detailed descriptions of the derivations of the tilted ring models and rotation curves of the individual galaxies. For 
each galaxy we show the Hermite velocity field, the resulting tilted-ring model velocity field, and the difference (residual) velocity 
field. Additionally we show overlays of the derived rotation curves on minor- and major-axis position-velocity (pV) diagrams, the 
rotation curves themselves, the radial distributions of PA and ;, and the positions of the dynamical center in comparison with the 
integrated H I maps and the SINGS 3.6 /im images'"*. We also briefly discuss the derivation of the various tilted ring parameters 
for each galaxy. 

The rotation curves, velocity fields and other supplementary information are shown in Figs. l68ll86l The following is a general 
description of the contents of these figures. 

Top row: Left: Velocity field derived from fitting Hermite polynomials to the natural-weighted data cube. See Sect. 13. II for 
a description. The systemic velocity is indicated by the thick contour The spacing AV between velocity contours is indicated 
in the figure. The approaching side can be identified by the light gray-scales and the dark contours. The receding side can be 
identified by the dark gray-scales and the white contours. The adopted dynamical center is indicated with a cross. The beam 
size is indicated by the ellipse enclosed by the rectangle in the bottom-left corner. Center: Model velocity field derived from 
the tilted-ring model. Gray-scales and contour levels are identical to those in the observed velocity field in the left panel. Right: 
Residual velocity field defined as the observed velocity field minus the model velocity field. The gray-scale range runs from -20 
kms"' (white) to +20 kms"' (black). Contour levels are -10,-20,-30, ... kms"' (black) and +10,+20,+30, ... kms"' (white). 

Centre row: Left: Position-velocity diagram taken along the average PA of the major axis as listed i n Tabled This PA is also 
indicated in the top-left of the panel. The thickness of the slice equals one pixel (typically 1.5", cf. IWalter et aDl2008h in the 
corresponding cube. Contours start at +2a in steps of 4ct (full contours), and -2a in steps of -4cr (dotted contours). The systemic 
velocity and position of the center are indicated by dashed lines. Over-plotted is the rotation curve projected onto the average 
major axis using the derived radial variations of PA and /. The spatial and velocity resolutions are indicated by the cross enclosed 
by the rectangle in the bottom-left corner. Right: Position-velocity diagram taken along the average minor axis. Contours and 
symbols are as in the major axis diagram. Over-plotted is again the rotation curve but projected onto the average minor axis using 
the derived radial variations of PA and /. 

Bottom-left panel: Top: The 3.6 /im IRAC Spitzer image. For all galaxies the same logarithmic intensity scale was used, 
running from log(//(MJy ster"')) = -2 (white) to log(//(MJyster~')) = +1 (black). The dynamical center is indicated by a cross. 
Bottom: Integrated H I map. The dynamical center is indicated by a cross. The maximum column density level displayed (black) 
is 2 • 10"' cm"^. The contour indicates the 3cr level. This level was computed as follows: for the standard THINGS cubes 
Ctot = VNaomn, where (Ttot is the noise in the integrated H I map, cTchan the noise in one channel, and the number of channels 
contributing to a pixel. 

Bottom-right panel: Top: The rotation curve corresponding to the tilted-ring model is represented by the black filled circles. 
The error-bars correspond to the dispersion of the velocity values found along the corresponding ring. This rotation curve was 
derived with all parameters fixed to their final values (as indicated by the thick black curves in the PA and inclination panels 
below). The full drawn gray line shows the corresponding rotation curve of the approaching side derived using these adopted 
distributions of inclination and PA. The dashed gray line shows the equivalent rotation curve of the receding side. The dotted 
black curve indicates the rotation curve of the entire disk derived with PA and inclination as free parameters. The uncertainties 
in the rotation velocity are defined as the quadratic addition of the dispersion in velocities found along each ring and one quarter 
of the difference between the approaching and receding sides velocity values. See Sect. 3.5 for a full discussion. Center and 
bottom: Inclination and PA values used in the tilted-ring models. The open circles show the distribution of inclination and PA 
when left as free parameters. These values result in the dotted black rotation curve described above. The crosses show same 
for the approaching side, the gray filled circles for the receding side. The thick black lines indicate the distributions for PA and 
inclination that were ultimately adopted to derive the rotation curves in the panel above. 

We now discuss details of the derivation of tilted-ring models of individual galaxies. 

NGC 925: The results for NGC 925 are shown in Fig. |68] NGC 925 is classified as a late-type barred spiral. A small weak 
central bar-shaped component is indeed visible in the IRAC 3.6 /im image. The outer parts are dominated by a two-armed spiral, 
with the southern arm also prominently visible in H I. With the dynamical center fixed, we find that Vsys varies between ^ 555 
kms"' in the center to ~ 545 kms"' a.t R ^ 120". At larger radii Vsys settles at a constant value of 546.6 ±0.6 kms"'. This 
value is very close to the average value determined over the entire radial range (546.3 ± 3.9 kms"'), and we adopt the latter 
value as the systemic velocity. This value is slightly lower than the central value of the global profile of 552.5 (551.5) kms"' as 
measured at the 20 (50) percent level. The PA distribution is continuous and well-defined. The inclination shows some scatter, but 
a global trend from / ~ 75° in the inner parts to / ^ 60° in the outer parts is clearly visible, indicating a slight warp in NGC 925. 
The rotation curve derived using the PA and inclination models is almost identical to the one with PA and ; as free parameters. 
Despite the asymmetric appearance of the H I disk, the approaching and receding sides are fairly symmetric in terms of their 
rotation velocity. 

NGC 2366: The data and analysis of NGC 2366 are given in lOh et aP (l2008l) . For completeness, we summarize and the data 
and analysis in Fig. |69l but refer refer to their paper for a complete discussion. 

NGC 2403: NGC 2403 is a late-type Sc spiral and member of the M81 group. We present the data in Fig.|70l The IRAC 3.6 
^m image shows multiple spiral arms, some of which are also evident in the integrated H I map. The velocity field of NGC 2403 
is very regular as shown by the perpendicular major and minor axes and the close agreement between approaching and receding 

Note that NGC 2366, NGC 2903 and DDO 154 are not part of the SINGS sample. For these galaxies we have retrieved data from the Spitzer archive. 
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sides. The small wiggles in the inclination distribution are caused by streaming motions induced by the spiral arms, also visible 
as residuals in the residual velocity field. After fixing the position of the dynamical center, we find the systemic velocity to be 
well-behaved, showing no large-scale trend with radius. The average value is Vsy, = 132.8 ±1.6 kms"', which is the value we 
adopt here. This value is close to the central velocity of 133.1 (134.8) kms"' derived from the global profile at the 20 (50)percent 
level. The inclination and position angle are well-behaved. A comparison of the curve with PA and / as free parameters with 
the curve which has these parameters fixed using the model distributions shows no appreciable differences. Note the remarkable 
symmetry between approaching and receding sides in rotation velocity, inclination and position angle between 450" ^R'^ 850". 

NGC 2841: NGC 2841 is an early-type (Sb) spiral which is dominated by a prominent central bulge, as well as a central hole in 
the H I distribution, as shown in Fig.|7Tl The tightly-wound spiral structure is also visible in the inner part of the H I distribution; 
in the outer parts the disk is dominated by two arms which seem to bend away out of the plane of the inner disk. The velocity field 
is regular in the inner parts but in the outer regions shows clear evidence for the presence of a warp. The major-axis slice shows 
that at the outer edge of the central H I hole the rotation has already almost reached its maximum value. After fixing the central 
position, we find a systemic velocity that is constant within R < 350", but starts to deviate at larger radii. Due to the obvious 
presence of the warp at these larger radii we adopt the average value within R = 350" as our best value: Vsys = 633 .7 ± 1 .8 km s~' . 
This agrees well with the central velocity derived from the global profile at the 20 (50) percent level of 635.2 (636.2) kms"'. 
With the systemic velocity and dynamical center fixed, we find well-defined distributions of PA and ; with radius. Despite the 
warp, NGC 2841 is very symmetrical with no obvious differences in PA and / for the approaching and receding sides. Our best 
model for the inclination shows an approximately constant value of ; ^ 70° for R < 200", with a gradual but systematic increase 
to / ~ 78° in the outer parts. The rotation curve with PA and / as free parameter is identical (within the uncertainties) to the one 
derived with PA and / fixed. 

NGC 2903: As can be seen in Fig. |72] the appearance of NGC 2903 is dominated by tightly-wound spiral arms and a bar as 
visible in the 3.6 fim image. The velocity field and integrated H I map can be divided into two distinct regions: in the inner parts 
we find an H I ring, corresponding to the bar and inner spiral arms as seen in the IRAC image. At the outer radii we see a lower 
surface density disk which is dominated by two broad H I spiral arms. These two regions are clearly separated around R ^ 240" 
by minima in the H I column density. The kinks in the velocity contours near the outer spiral arms clearly indicate the presence of 
non-circular motions there. The very inner parts of the disk also show a strong twist in the PA of the kinematical major axis, very 
likely due to the effect of the bar and associated streaming motions. Their presence is also implied by the motions that can be 
seen in the minor axis profile and the increased width of the profiles in the inner parts (visible in the major axis position-velocity 
diagram). With the position of the dynamical center fixed, the systemic velocity still shows large changes in the inner parts, most 
likely due to the effects of the bar We fixed Vsys to the average value between 200" <R < 400" and find Vsys = 555 .6 ± 1 .3 km s~' . 
This agrees well with the central velocity from the global profile derived at the 20 (50) percent level of 556.6 (555.7) kms"'. 
The inclination and PA for R > 100" are straightforward to model, but at smaller radii both show large variations. The sudden 
change in inclination at R ^ 110" corresponds to the outer radius of the bright, compact spiral arms associated with the central 
bar, whereas the sudden change in PA at ~' 60" occurs where the bright tips of the straight central bar are found. As both of 
these changes are associated with non-rotational motions, we have for our tilted-ring models simply extrapolated the innermost 
unaffected values inwards. 

NGC 2976: NGC 2976 is classified as an Sc galaxy. Its analysis is illustrated in Fig. [73] The IRAC 3.6 fim image shows 
no sign of spiral arms or bar component. The stellar component of the galaxy is characterized by two regions of enhanced star 
formation at either end of the disk, which correspond to the two density enhancements seen in H I. Having fixed the dynamical 
center, we do not detect systematic trends of Vsys with radius except for a slightly larger scatter in the innermost few points. We 
therefore fix it to the average value for/? > 20" and find Vsys = 1.1 ± 1.3 kms"', which compares well with the central value of 
2.6 (4.2) kms"' derived from the global profile at the 20 (50) percent level. The trend of PA with radius is well defined, and can 
can be approximated by a constant PA for R > 70", with a modest increase within that radius. The inclination shows increased 
scatter within R < 70", and we fix it to the average value derived for 80" < R < 1 10". The final rotation curve agrees well with 
the curve derived with PA and / free. The residual field also shows no large residuals, indicating that circular motions describe 
the dynamics of this galaxy well. 

NGC 3031: The analysis of NGC 3031 in illustrated in Fig.|74l The IRAC 3.6 ^im image as well as the H I total intensity map 
show two well-defined spiral arms. The H I observations used here consist of a mosaic of two separate pointings. Because of this, 
it was necessary to use a primary-beam corrected data cube. These corrections cause an apparent increase of the noise towards 
the edges of the field, apparent in the major- and minor-axis slices. Careful blanking of these regions ensured that they did not 
affect the analysis presented here. Also note that due to the very large velocity width of the gas in the M81 system it was not 
possible to define line-free channels in order to subtract the continuum (seelWalter et al. 2008). The total intensity H I map, as 
well as the position velocity slices therefore still contain the radio continuum emission. The central radio continuum source can 
be seen as the vertical high-intensity feature in the position- velocity slices. An inspection of the velocity field shows that M81 
has clear non-circular motions associated with the prominent spiral arms. The residual velocity field clearly shows the effects of 
these streaming motions. It is also clear that beyond R > 800" the H I gas is no longer in regular rotation around M81. At these 
large radii the motion of the gas is starting to become dominated by the tidal interaction processes within the group. For this 
reason we do not attempt to derive the rotation curve beyond that radius. More evidence for these tidal effects is evident in the 
minor axis slice where we see a sudden and large deviation of the gas from the systemic velocity. This corresponds to the strong 
twist in the systemic velocity contour seen in the eastern part of the outer velocity field. After fixing the position of the dynamical 
center, we determine the systemic velocity by first fixing ; and PA to indicative values of 60° and 330°, respectively. We find Vsys 
to be constant within R < 350" and use the average value of -39.4 ±2.8 kms"' . This agrees with the central value of the global 
Hi profile of -39.4 (-41.3) kms"' at the 20 (50) percent level. After fixing Vsys, and re-running ROTCUR we find well-defined 
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distributions of / and PA. The radius of ^ 400" at which the sudden change in inclination occurs corresponds to the outer edge of 
the innermost H I "ring" visible in the integrated H I map. From the PA distribution, it is clear that the approaching and receding 
sides are not entirely symmetric. The rotation curve itself, by contrast, is. The curve shows a very pronounced drop in the outer 
parts; however given the position of M81 within an interacting system, this very likely carries no implication for the distribution 
of the dark matter. 

NGC 3198: The analysis of NGC 3031 in illustrated in Fig.|75] The IRAC 3.6 //m image shows two well-defined strong spiral 
arms, with several less prominent arms branching off from the two main ones. The velocity field is regular with indications of a 
very modest warp. Small streaming motions due to the two dominant spiral arms show up as small kinks in the velocity contours. 
These are also visible in the residual velocity field. With the position of the center fixed, we find a well-behaved systemic velocity 
with an average value Vsy, =660.7 ±2.6 kms~', which compares favorably with the central velocity of 661.2 (661.7) kms"' as 
derived from the global profile at the 20 (50) percent level. The PA is well-defined for all radii; for R > 300", the distribution 
of i is also well-defined, and shows a small but systematic increase. For R < 300" the inclination shows some variations due to 
motions along the spiral arms (cf. the low-amplitude pattern visible in the residual velocity field). We therefore adopt an almost 
constant inclination for R < 300". Note that the rotation curve does not depend on these assumptions: the curve with ; and PA as 
free parameters equals the final curve almost point by po int. The approaching and receding sides are very symmetrical. 

IC 2574: The data and analysis of IC 2574 are given in lOh et al. (2008). For completeness, we show a summary of their results 
in Fig.|76] but refer to their paper for a complete discussion. 

NGC 3521: NGC 3521 is a disk galaxy showing a flocculent spiral structure. The inner Hi disk seems to be in regular 
rotati on, but the outer disk shows a much different behavior This is most clearly visible in the channel maps (see IWalter et alJ 
l2008h . but can also be seen in the major-axis position-velocity slice, as shown in Fig.|77] The central velocity derived from the 
symmetrical global profile is 798.2 (798.6) kms"' at the 20 (50) percent level. After fixing the position of the dynamical center, 
the systemic velocity shows a systematic trend with radius and the largest radial variation observed within the THINGS sample. 
It is constant at a value of ^ 800 kms"' for R < 250" (8 kpc) and R > 450" (15 kpc), but drops to a value of ^ 111 kms"' at 
intermediate radii. Due to the good agreement between Vsys values found in the inner and outer parts, we fix Vsys to the average 
value for R < 250" and R > 450" and get Vsys = 803.5 ± 4.5 kms"' . This choice for Vsys leads to a rotation curve which shows a 
large difference between approaching and receding sides at the intermediate radii. 

NGC 3621: NGC 3621 is a late-type spiral with an Hi d istribution that is c haracterized by a regularly rotating disk and a 
warp in the line of sight (see the channel maps presented in I Walter et al.ll2008h . This feature is also visible in the major-axis 
position-velocity diagram, as shown on Fig.|78] To derive the velocity field of the disk component only, we blanked all emission 
associated with the warp prior to deriving the velocity field. This velocity field and the associated integrated H I map presented 
here, therefore only show the disk component. For the full velocity field and H I map see i Walter et al.i (j2008i) . After fixing the 
position of the center, we find Vsy, to be on average constant with radius. Small variations with a maximum amplitude of 5 
kms"' are visible, but show no systematic trends with radius. We adopt the average value over the entire radial range as our best 
value: K^ys = 728.5 ±2.7 kms"' . This is consistent with the central value of 730. 1 (728.8) kms"' as derived from the global H I 
profile at the 20 (50) percent level. The PA is well-defined and only varies slightly with radius. The inclination shows significant 
variations, but in general decreases with radius. 

NGC 3627: NGC 3627, or M66, is a barred, interacting spiral galaxy that is part of the Leo Triplet. Our analysis of this 
galaxy is illustrated in Fig.|79] The IRAC 3.6 /im image shows a pronounced bar, with a two-armed, asymmetric spiral structure. 
The western arm is also clearly visible in the integrated Hi map. The eastern arm is confused with other high-column density 
features. The pronounced spiral arms and asymmetry of NGC 3627 do not make this a prime candidate for detailed studies that 
assume azimuthal symmetry. In the central parts the H I column density lies below the 3(7 level that we applied during the profile 
fitting, hence the absence of velocity information there. After fixing the position of the dynamical center, we fix PA and ; to 
indicative values of 175° and 60°, respectively, and find Vsys values that vary slightly within R ^ 100", but are constant outside 
that radius. We use the mean value for R > 100" and find Vsys = 708.2 ±1.1 kms"'. This differs somewhat from the 20 (50) 
percent level central value derived from the global Hi profile of 717.3 (720.3) kms"', but given the asymmetry of this galaxy 
this should not come as a surprise. Trends of / and PA with radius are fairly straightforward to determine with both parameters 
showing modest variations. The approaching and receding sides show different radial trends, as might have been expected from 
the obvious asymmetry visible in the IRAC image. Nevertheless, the respective rotation velocities converge in the outer parts, 
leading to a well-determined flat outer rotation curve. 

NGC 4736: NGC 4736 is dominated by an inner high surface brightness disk in the IRAC 3.6 fim image as shown in Fig.lSOl 
Prominent arms are visible in the integrated H I map. After fixing the position of the dynamical center, the systemic velocity 
shows no systematic trends with radius and we adopt the average value Vsys = 306.7 ± 3.7 kms"'. This corresponds well with 
the values derived from the global profile at the 20 (50) percent level of 307.6 (306.6) kms"'. The PA is very well defined and 
constructing a model is straight-forward. The inclination shows a well-defined trend from ~ 45° in the outer parts to ~ 35° 
at R ^ 80". At this radius the slope of the inclination steepens, and the value appears to decrease to ~ 20° in the center We 
consider this break in the inclination slope not physical; the inner region of NGC 4736 is known to be a region of non-circular 
motions and the radius w here the break occurs coin cides with the radius of the inner star formation ring (the "expanding Hll 
regions ring" described in iMunoz-Tuiion et al.ll2004l) . The minor-axis position-velocity diagram also shows clear evidence for 
strong non-circular motions in the inner arc-minute. For our inclination model we simply extrapolate inward the trend found 
for R > 80". Despite the disturbed distribution of the H I, the approaching and receding sides of the final rotation curve are 
remarkably similar. The velocity systematically decreases with radius: it peaks at V = 21 1 kms"', and declines to V = 1 15 kms"' 
in the outerr nost parts. As with NGC 3627, it is clear from the morphology that the assumption of azimuthal symmetry is not 
appropriate. iTrachternach et all ( l2008h show that NGC 4736 has the second highest non-circular motions in the current sample 
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(after NGC 3627) and any interpretation of its rotation curve should therefore be treated with care. 

DDO 154: Our analysis of DDO 154 is illustrated in Fig.[8Tl DDO 154 has no central radio continuum source or central light 
concentration to help pinpoint the dynamical ce nter. The position of the dynamical center was therefore derived purely using the 
velocity field (see also lTrachternach et alj|2008l) . The systemic velocity shows a stable behavior. Using the average value within 
R = 260" we find Vsys = 375 .9 ± 1 .4 km s"' . This corresponds very well with the central values derived from the global H I profile 
at the 20 (50) percent level of 375.5 (375.4) kms"'. The PA is straightforward to model. The inclination shows some variation in 
the inner parts, mainly due to the small number of independent pixels per ring; here we simply extrapolated the model inwards 
from the smallest radius where / was well-behaved. The outer parts show an apparent steep drop in inclination. Most of the gas 
contributing to these points is however found in the outermost flocculent parts of the warp, and given the small filling factor of 
the rings there, it is not clear how realistic this steep drop in inclination is. We therefore extrapolated our model outwards from 
the largest radii with well-determined inclinations. 

NGC 4826: As can be seen in Fig. [82l the outer H I disk of NGC 4826 has a very low column density - one of the lowest in 
the entire THINGS sample. Only a small fraction of the disk has profiles with a peak value > 3(t. Nevertheless, this small filling 
factor, combined with the very regular kinematics in the outer disk is still sufficient to derive a well-defined outer rotation curve. 
Small insets in Fig. [82] show the region of the counter-rotating disk, as well as the corresponding model and residual velocity 
fields in more detail. After fixing the position of the dynamical center, we find no systematic trends of the systemic velocity with 
radius, and adopt the average value: Vsys = 407.4±7.0 kms~'. This value corresponds very well with the central value of 407.9 
kms"' as derived from the global H I profile for both the 20 and 50 percent level. The inclination and PA show clear trends with 
radius. A good description for the inclination is a constant value of ~ 55° in the inner parts, followed by break and a sudden 
increase to ~ 70° atR^ 180", in turn followed by a gradual decline in the outer parts to a value of ~ 63°. 

NGC 5055: As illustrated in Fig. [83] the major-axis position-velocity diagram shows evidence for the presence of gas at 
velocities lower than the local rotation velocity. The structure of NGC 5055 as seen in the 3.6 ^m image is regular, with signs of 
well-defined flocculent spiral arms. After fixing the position of the dynamical center we find that the systemic velocity shows a 
systematic variation, from a constant value of ^ 496 kms"' within R < 200" to a value of ^ 510 kms"' in the outer parts. We 
make the assumption that the gas in the inner parts is likely to be more tightly bound to the system than the tenuous outer arms, 
and adopt the average value between 50" < R < 180". We find Vsys = 496.8 ±0.7 kms"'. The central velocity as derived from the 
global profile at the 20 (50) percent level is 499.3 (497.4) kms"' . The PA and ; values are well-defined for the inner H I disk. The 
distributions for the approaching and receding sides of the outer warped part differ significantly, leading to different approaching 
and receding rotation curves for R > 250". 

NGC 6946: Our analysis of NGC 6946 is illustrated in Fig. [84] For the position of the dynamical center we adopt the 
position of the central radio continuum source. The systemic velocity is fairly constant with radius. We adopt the average value 
within 500" of Vsys = 43.7 ± 3.3 kms"'. This agrees well with the central value of 42.0 (45.1) kms"' derived from the 20 (50) 
percent level of the global H I profile. The PA is well-behaved and straight-forward to describe. The inclination, when left free, 
shows a large scatter, showing the increased difficulty of determining the inclination using only dynamical information at low 
inclination values. Nevertheless, the large number of independent data points allow us to accurately model it as a linear radial 
dependence. The assumption we make is that the scatter around this average value is caused by small non-circular motions that 
have a disproportionally large influence due to the low inclination. Due to this inclination scatter we do not consider here the 
rotation curve derived using PA and / as free parameters. 

NGC 7331: We illustrate our analysis of NGC 7331 in Fig. [85] The velocity field, though regular on large scales, shows many 
wiggles in the velocity contours, often coinciding with the positions of the spiral arms. After fixing the position of the dynamical 
center, we find that the systemic velocity is well-defined over the range 80" < R < 280". Between 80" and 200" it varies 
systematically from ~ 830 kms"' to ~ 805 kms"', settling on a constant value of 818.3 ±0.9 kms"' between 200" < R < 280". 
We have chosen the latter value as our Vsys value, comparing favorably with the central velocity from the global profile of 
815. 7 (816.5) kms"' as measured at the 20 (50) perc ent level. The varying Vsys value could indicate a kinematic lopsidedness 
("e.g.. ISchoenmakers et al.iri997t ISchoenmakerill999l) . A harmonic analysis of the velocity field to quantify this is presented in 
iTrachternach et al.l (l2008h . In the outer parts the PA values for approaching and receding sides differ distinctly. At R ^ 150" 
the inclinations are different as well, most likely due to streaming motions along the spiral arms. These show up clearly in 
the residual velocity field. Also note the large difference between the rotation curves of approaching and receding sides within 
R < 200". The typical difference is ~ 25 kms"', and shows why V^ys is difficult to constrain over this radial range. The major 
axis position-velocity diagram shows some Hi outlining the rising part of the rotat ion curve, but the emission is too faint to 
fit a rotation curve. NGC 7331 was previously observed in Hi bv iBegemanI ( 19871) . who already noted many of the features 
we have describ ed above. The dynamical center as determined by Be gemanl (Il987h . coincides with ours to within the errors. 
lBegemanl(ll987h also notes the same varying Vsys in the inner parts. He determines Vsys using only points at > 210" and finds 
a value 820.1 ± 2.7 kms"'. For the P A a trend very similar to ours was found. The trend we find in the inclination distribution 
is also present in the lBegemanl(ll987h data, but with the much smaller number of resolution elements in that data set, a constant 
inclination model was chosen there. 

NGC 7793: As shown in Fig. [86] the velocity field of NGC 7793 is well-defined with indications of a small-amplitude warp. 
The major-axis slice shows little evidence for extensive non-circular motions. The residuals are small, and prominent in only a 
few small areas in the outermost part of the disk. After fixing the position of the dynamical center we find the systemic velocity 
to be well-defined: we adopt the average value between 40" < R < 250" and find Vsys = 226.2 ± 1 .2 kms"', agreeing well with 
a central velocity of 227.2 (226.6) kms"' as derived from the global profile at the 20 (50) percent level. The PA and ; values 
are straight-forward to model: the PA is well defined; the inclination shows a modest increase in the inner parts, followed by a 
gradual decline in the outer parts. The curve as derived with ; and PA as free parameters, is within the uncertainties equal to the 
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one with / and PA fixed to their model values. 
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Fig. 68.— 



Summary panel for NGC 925. See the Appendix for more information. 
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Fig. 69.— 



Summary panel for NGC 2366. See the Appendix for more information. 
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Fig. 70.— 



Summary panel for NGC 2403. See the Appendix for more information. 
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Fig. 71. — Summary panel for NGC 2841. See the Appendix for more information. 
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Fig. 72.— 



Summary panel for NGC 2903. See the Appendix for more information. 
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NGC2976 





Fig. 73.— 



Summary panel for NGC 2976. See the Appendix for more information. 
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Fig. 74. — Summary panel for NGC 3031. See the Appendix for more information. 



THINGS rotation curves 



93 



NGC3198 




Fig. 75. — Summary panel for NGC 3198. See the Appendix for more information. 
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Fig. 76. — Summary panel for IC2574. See the Appendix for more information. 
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Fig. 77.— 



Summary panel for NGC 3521. See the Appendix for more information. 
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NGC3621 





Fig. 78 



. — Summary panel for NGC 3621. See the Appendix for more information. 
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NGC3627 




Fig. 79. — Summary panel for NGC 3627. See the Appendix for more information. 
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Fig. 80.— 



Summary panel for NGC 4736. See the Appendix for more information. 
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Fig. 8 1 . — Summary panel for DDO 154. See the Appendix for more information. 
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Fig. 82. — Summary panel for NGC 4826. See the Appendix for more information. 
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Fig. 83. — Summary panel for NGC 5055. See the Appendix for more information. 
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Fig. 84.— 



Summary panel for NGC 6946. See the Appendix for more information. 
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Fig. 85.— 



Summary panel for NGC 7331. See the Appendix for more information. 
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NGC7793 




Fig. 86. — Summary panel for NGC 7793. See the Appendix for more information. 
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